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ABSTRACT 

The distribution of inventory demand over procurement 
leadtime is modeled using mixed probability distributions 
that explicitly account for the high incidence of zero 
demands observed in Inventory Control Point Demand History 
Files. Analysis was limited to the right-hand tail area of 
the distribution on the assumption that that area is of crit- 
ical importance in determining the effectiveness of an inven- 
tory system. Probabilistic models studied were: 1) 

Bernoulli-exponential/ 2 ) Bernoulli-lognormal, and 3) 
Bernoulli-logistic, These compound distributions were com- 
pared to several standard distributions including the Pois- 
son, negative binomial and normal distributions using a 
resampling procedure appropriate in cases such as this where 
a paucity of data exists. Fits obtained from the 75th to 
95th percentiles indicated the mixed distributions may be 
superior as a class to the standard distributions for high- 
demand items. 
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I. INTIO^CTIDN 



Inventory system performance^ in a situation with random 
demand depends upon reorder points computations. Reorder 
point computations depend upon the probabilistic model chosen 
to represent the inventory system in that the reorder point 
is composed of two parts: 1) The expected demand during a 

procurement leadtime and 2} The safety level^ determined from 
the probabilistic model of demand over leadtime. Probabilis- 
tic models are generally only utilized to represent demand 
over procurement leadtime since That is the only time period 
when stockouts^ potentially occur. The chosen distribution 
is then a conditional distribution given procurement lead- 
time. The Navy has historically used three distributions of 
demand for models utilized at the Inventory Control Point 
(ICP) echelon of the Navy Supply System. The current proba- 
bilistic models and the average annual demand that is used to 
determine which model to utilize are displayed in Table l.s 



1 See Appendix A for definition of '’System Performance". 

2See Appendix A for definition of "Reorder Point". 

3See Appendix A for definition of "Safety Level". 

♦See Appendix A for definition of "Stockouts". 

sThe Aviation Supply Office (ASO) has recently changed to 
the use of the normal distribution as the model of choice for 
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Table 1: Current Inventory Control Point 

Probabilistic Demand Over Leadtime Models 



D istr ibut i pn 



Average Annual 
2 § 1 § nd Ra nge 



Poisson 



0-1 Low 



Negative Binomial 1 - 20 Medium 



Normal 



> 20 High 



If the actual quarterly demand records for an ICP are 
examined, one striking characteristic that is present in 
nearly all records is the high incidence of zero observa- 
tions. This is not surprising when one considers the ech- 
eloning of the supply system and the facr that users only 
place a demand on the wholesale system as their reorder 
points are reached. Since most activities stock at least one 
quarter's expected demand for an item, reorders are expected 
at most four times a year. By the natural phase differences 
of each activity's reorder actions, the ICP may experience 
zero demands during a given quarter for even high-demand 
items, and will surely often experience quarters of zero 
demands for medium- and low-demand items. A previous study 
of Pacific Fleet Combat Stores Ships' demand data [Ref. 1] 
found that the frequency of zero observations in those demand 



all demand categories. Ships Parts Control Center (SPCC) 
still uses all the distributions in Table 1. 
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records appeared to have a low correlation with the overall 
demand level (i.e. high, medium or low).* 

A desirable model for an inventory system should compute 
reorder levels that accurately correspond to a specified 
stockout risk based on the level of actual or anticipated 
demand for each item. If such a model were available, the 
inventory system could avoid the added expense of either hav- 
ing too much or too little stock. One problem with the dis- 
tributions in Table 1 is that none can accounr for the 
probability mass at zero; therefore, none can accurately com- 
pute reorder levels. 

This study explores the use of compound probability mod- 
els and their potential to more accurately achieve the goals 
presented above by explicitly accounting for, the probability 
mass at zero. The results obtained should be considered as a 
first step in the exploration of a class of models which have 
heretofore received little attention in the context of model- 
ing inventory demand. 



*See Section II.B.1 for a more detailed discussion of 
this study. 
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II. 



BACKG^UND 



A. CURRENT MODEL SELECTION TECHNIQUES 

In the past, the usual approach of choosing the probabi- 
listic model for lead time demand was based on matching the 
empirical cumulative distribution function determined from 
the actual demand to various theoretical cumulative distribu 
tion functions. The shapes of the two curves are compared 
using the Kolmogoro v-Smirno v Statistical Goodness of Fit 
test. 7 The Chi-Squared test has historically not been used 
because it is not very powerful when the number of observa- 
tions is small {supply demand data is retained for at most 
twelve quarters and usually only eight quarters at the ICP 
level). The Kolmogorov-Smi rnov test provides only a relativ 
measure of goodness of fit, but it is usable with very small 
sample sizes. The goal then is to select the model with the 
best agreement between the empirical and theoretical cumula- 
tive distribution function curves. 

The Navy Fleet Material Support Office (FMSO) has uti- 
lized the above Kolmogorov- Smirnov test to evaluate several 
theoretical distributions. Most recently [Ref. 2] FMSO 

^A good description of the Kolmogorov-Smirno v one and 
two-sample tests may be found in Siegel, S., Hon Parametric 
Statistics , p. 47-52 and 127-136, McGraw Hill, 1956. 
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evaluated seven distr ibu tio nss at the 90*^ confidence level 
using demand data from NAS Brunswick, The use of demand data 
from a stock point vice an ICP assumes the underlying demand 
behavior at the stock point and ICP are similar. Though this 
was never validated directly, the NAS Brunswick demand data 
were deemed to have "unique qualities better suited for 
demand analysis" as follows; 

1. Twelve quarters of demand and demand frequency data 
were available per item, 50^ more than available in ICP 
files, 

2. New item identifiers that permit the selection of only 
steady state demand items were available, 

3. Reliability of the data was established from prior FMSO 
studies . 

The FMSO analysis showed that Navy demand patterns are 
very poorly modeled by any of the seven distributions stud- 
ied. No new models were proposed and no attributes of the 
demand patterns were indicated as the most significant factor 
in the failure of the standard distributions to model Navy 
demand patterns. 

There are several problems winh the way the Kolmogorov- 
Smirncv test was used in the FMSO study: 



®Distr ibutions tested were: Normal, Negative-Binomial, 
Poisson, Logistic, LaPlace, Gamma and Uniform. 
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1. The procedure is not strictly applicable in the case of 
discrete distributions such as the Poisson and negative 
binomial. 

2. The test assumes that the parameters of the distribu- 
tion being tested are completely specified in advance, 
i.e. not estimated from the data. (Lilliefors [Ref. 3] 
developed alternate tables for use in testing the expo- 
nential distribution when parameters are estimated from 
the data. This procedure was used in Reference 1.) 

3. The Kolmogorov- Smirno V test evaluates goodness-of-f it 
over the entire range of the distribution. For inven- 
tory problems, the region of greatest interest is the 
right-hand tail of the distribution. 



3. PREVIOUS RESEARCH 

Research into distributions which have properties that 
more closely model actual demand patterns has been very lim- 
ited and, when conducted, often did not explore anything but 
the usual standard distributions. Two projects though stand 
out because of their innovative use of non-standard distribu- 
tions and the good fit achieved. 

1 . Ber noulli-e xp onent i al D i strib utio n 

In Reference 1, Pacific Fleet Combat Stores Ship 
demand was modeled using the mixed Bernoulli-exponential 
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distribution. This distribution was chosen because it was 
observed that nearly all line items had an unusually high 
occurence of demands of zero. Analysis of the data indicated 
that with high probability, the number of non-zero observa- 
tions for a given item was unrelated to the average value of 
those non-zero observations. The conclusion was that the 
demand process could be thought of as two independent subproc- 
esses with one process determining whether a demand would 
occur or not and the other determining the quantity of the 
demand, given the demand did occur. The former process was 
modeled as a Bernoulli process with parameter p, p being the 
probability that a demand did occur. The latter process was 
modeled as a continuous exponential process based on explora- 
tory data analysis. The resulting probability distribution 
function for the number of units demanded in a leadtime for 
the Bernoulli -exponential distribution is given by: 



1 - p , X = 0 

h(x) = ^ 



p Xex p (- Xx) , X > 0 , 



( 1 ) 



where : 

VX = the expecred value of demand, given that the 
demand is greater than zero, 
p = probability that demand does occur. 
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Th9 

by; 



complementary cumulative distribution function is given 



H <x) = p exp (-Ax) , X > 0, (2) 

where H (x) is the probability that demand will not exceed x; 
this is equivalent to the complementary cumulative distribu- 
tion function used by Hadley and Whitin [Ref. 4]. 

The hypothesis that the demand of Combat Stores Ships 
could be modeled according to the above probability distribu- 
tion was tested using the X olmogorov- Smirnov test described 
earlier plus tests on the theoretical risk and the observed 
risk for values of risk in the upper right hand tail area of 
the distribution. The results of these tests on five inde- 
pendent samples from the demand data available provided 
strong evidence that the Bernoulli-exponential distribution 
describes the demand for any given stock item very well. 

2 . Loqar it hmic-Poisso n -Gam ma D is t r ibution 

A model for the distribution of demand during pro- 
curement leadtime [Ref. 5] was derived under the following 
assumptions ; 

a. Requisitions occur according to a stationary Poisson 
Process, 

b. Requisition sizes follow a logarithmic distribution. 
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c. Leadtime is a random variable with the gamma 
distribution. 

The resulting probability function for the number of units 
demanded in a leadtime is: 

h(x) 

and 
h(0) 
where ; 




cx,S = parameters of the Gamma distribution defining 
leadtime, 

0 = parameter of the logarithmic distribution 
solved for by interval bisection from 



S (X) = -0 / (1 - 0) la (1 - 0) , 
X = requisition arrival rate, 
c = - X / ln(1 - 0) 



^x,k 



Stirling numbers of the first kind 
from the recursion 

^x,k“ ^x-i,k-l * ~ ^)^x-l,k 



computed 



for k = 1, 2, ..., X and x = 1, 2, ... 

with Sv r, =3 for all x 
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This distribution is called the Logarithmic-Poisson-Ganiaa 
(LPG) . It is defined from four parameters a, 6, 0, and X 

(0 and X determine c) . 



For more efficient computations, a recursion equation 
for integer a was developed. 

X 



and 0, c, a, 6 and X defined as in (3). 

Computing reorder points from fractiles of the LPG 
distribution, even using the simplified version of (4), is 
generally too time consuming for many real applications. As 
a result, an approximation was developed and tested that used 
a scaled version of the Poisscn distribution to approximate 
the negative binomial distribution. The resulting approxima- 
tion is : 




(4) 



k=i 



where 



0 




p(Z(t) = kx} 




for X = 0, 1, 2, 



whe re 



Z(t) = number of units demanded in time t 
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M = C 0 



ct, g, c and 0 defined as in (3). 

The authors conducted subjective tests of goodness of 
fit to their LPG model and its derivatives using Air Force 
consumable items and concluded, using as an example one item, 
that there was a "very close agreement between the observed 
and the predicted cumulative distribution functions for this 
item. " 
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III. PROPOSED MODELS 

A. COMPOUND BERNOULLI MODELS IN GENERAL 

The results of the previous research indicate that a 
potential exists for improving inventory performance by the 
use of slightly more complex models. The focus of this study 
is on the compound Bernoulli models, leaving further analysis 
of the LPG model and its derivatives to others. The use of 
compound Bernoulli models has several advantages over stan- 
dard simple distributions: 

1. The model, by necessity, has an added parameter whose 
sole function is to estimate the probability of a non- 
zero demand. 

2. For simple models defined only over the positive real 
numbers, the addition of the Bernoulli parameter simply 
results in the scaling of the cumulative distribution, 
thus adding little additional complexity to the model. 

3. By explicitly accounting for zero demands, compound 
models account for more of the observed variance than 
do simple models. 

Compound models are not without disadvantages: 

1. Each parameter added to the model must be estimated for 
each individual stocSc item. For a large inventory 
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system such as the Navy's, this equates to considerable 
additional computational requirements. 

2. Additional parameters must be stored for future 
retrieval or else computed from stored data when 
required, thus either extending computational time or 
else requiring additional online storage space. 

The Bernoulli compound models all have similarly shaped 
cumulative distribution functions. Figure 1 is a sample dis- 
tribution function for the compound Bernoulli-lognormal dis- 
tribution with a Bernoulli factor of p equal to 0.6, a mean 
of 1.0 and a standard deviation of 0.5. Note that the dis- 
tribution has a mass of zero for negative observations, a 
mass of 1-p at zero and the usual cumulative distribution 
function shape for positive observations but of mass p vice 
mass one. 

B. aODEL DESCRIPTIONS 

Three compound Bernoulli models were formulated for eval- 
uation. Each model was derived from the base distribution 
after compounding with the Bernoulli process. An additional 
requirement was made that the inverse of the compound model 
must be computationally relatively simple. 

1 . Bernoulli-e xpon enti al D ist r ibuti on 

This distribution was derived in Reference 1 and the 
derivation will not be repeated here. Equations (1) and (2) 
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Compound Cumulative Distribution Function, 
Bernou lli-lognor mal Distribution 



are the probability density and complementary cumulative 
distribution functions respectively. The Bernoulli- 
exponential inverse complementary cumulative distribution 
function is computed as: 






0 



, H(x) > p 



X = < 



M[ln(p) -ln(H(x))], H(x) <p 



( 6 ) 



The parameters of this distribution are estimated from the 
sample data as follows: 
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p = number of non-zero demand observations in sample 
divided by sample size, 

V* = X / p, where x is the entire sample mean. 



or 

y 



X, where x is the sample mean of only the non-zero 
demand observations. 



2. Bernoulli -l ognorma l Dis t ribu tion 

The lognormal distribution was chosen as a candidate 
for testing because of several desirable properties: 

a. It is defined for all real numbers greater than zero, 

b. Its inverse is readily computed using inverse normal 
approximations that are well documented for accuracy, 

c. Alone, the lognormal distribution does not allow the 
case of zero demand, but when combined with the Ber- 
noulli distribution, the compound Bernoulli-lognormal 
distribution defines all the expected demand values 
greater than cr equal to zero. 

Dtilizing the Bernoulli parameter p as a scaling fac- 
tor for the lognormal density function, the compound 
Berno ulli- lognormal density function may be expressed as; 



h(x) 




X = 0 



X > 0 



(7) 



where : 
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p = probability of a Don-zero demand, 

M = mean of the natural log of all demands greater than 
zero , 

o = Standard deviation of the natural log of all demands 
greater than zero, 

= Standard normal density function. 

Integrating (7) results in the complementary cumulative dis- 
tribution function given by: 

P , X = 0 (8) 

H(X) = 

In (X) - w 

a 




X > 0 



where : 

<J> = Standard normal cumulative distribution function. 

The Bernoulli-lognormal inverse cumulative distribution func- 
tion is easily derived from (8) as: 



0 

X 



, H(x) > p (9) 



expfu + (1 - H(x)/p) ], H (X) < p 

The parameters are estimated from the sample data in the 
usual manner: . 

p = number of non-zero demand observations in sample 
divided by sample size. 
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y = X, vfhere x is the sample mean of the natural log of 



all non-zero demand observations 



o = Sf where s is the sample standard deviation of the 
natural log of all non-zero demand observations. An 
alternate method for estimating the mean and stan- 
dard deviation of the lognormal distribution is to 
use the method of moments. The procedure in effect, 
transforms the mean and standard deviation obtained 
from the untransformed sample data vice transforming 
the data first and then computing the sample mean 
and standard deviation. This method involves fewer 
logarithms and will run faster on a computer. The 
procedure is described in Appendix D. 

3. Ber noulli -l cqis r ic , Distribut ion 

The logistic distribution is a pseudo-normal distri- 
bution that is similar in shape to the normal, but is easier 
to handle mathematically. Its density function is defined as 
[Ref. 6]: 



while its cumulative distribution function is defined as: 




( 10 ) 



f (X) 



F(x) 
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where : 



U = mean of the entire sample, 
a = standard deviation of the entire sample. 

With the logistic distribution, it is no longer a simple mat- 
ter of scaling the distribution by the Bernoulli p parameter, 
since this distribution is defined over the entire real line 
and not just on the positive half. Consequently, the distri- 
bution requires a new constant of integration to replace the 
1/2 used in (11). Integrating (10) from 0 to ® provides the 
new constant as: 



which, when compounded with the Bernoulli distribution, gives 
the Bernoulli-logistic distribution density function as: 




oo 



( 12 ) 



/ 1 - p 



X = 0 (13) 



1 




h (X) = y p 




/ 



X > 0 



1 + ta nh 



, otherwise. 



where : 



p = mean of all non-zero demand observations. 
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a = standard deviation of all non-zero demand 
observations . 

After integrating (13), the compleaentary cumulative distri- 



bution function may be. written as; 



H(x) = p< 



1 - 



1 


- - \ 


1 + tanh/ 


> 

Try 


\ 








/ 7T(X - U )' 

for X > 0 



(14) 



The inverse cumulative distribution is: 






0 



, H(x) > ? (15) 






\/3c 



2d 



li + In^ 

^ H (x) 



V 



/ Try M - 1 y , H (X) < p 

1 > tanh( ) 

\2y/foJj 



The parameters are estimated from the sample data in the 
usual manner: 

p = number of non-zeco demand observations in sample 
divided by sample size, 

✓V _ ^ 

y = X, where x is the sample mean of all non-zero demand 
observations , 

a = s, where s is the sample standard deviation of all 
non-zero demand observations. 
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17. EVALUATION PROC BDURB 



A. THE DATA 

In order to test the validity of the proposed models, 
samples of actual demand were obtained from the Operations 
Analysis Department at FMSO. The data, accumulated from the 
demand history files of the Aviation Supply Office, was orig 
inally used as input data for the 5A (Aviation Afloat and 
Ashore Allowance Analyzer) [Ref. 7]. The data consists of 
1587 consumable 1R cog items and 2892 non-program-related^ 
repairable 2R cog items. The information for each item is 
contained in one master record and several subrecords. The 
master record contains identifying information on each irem 
such as the national stock number, replacement price, etc. 
Each subrecord contains up to forty-six demand records, each 
including the demand quantity and the day of the demand. A 
complete record layout is contained in Appendix E. To ease 
processing, the Julian Dates of the original demands were 
replaced by FMSO with a sequential date ranging from 1 to 
1500 representing approximately four years of available 
history. 



^See Appendix A for definition of "program related" 
ite ms . 
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For this application, additional screening was conducted 
to group the demands into thirty-day buckets, thus providing 
a demand time series of forty-eight observations. Further 
editing was required prior to analysis to remove the negative 
demands noted for many of the items. It was assumed that the 
negative observations were the result of cancellations of 
previous demands which were never filled. Therefore, when- 
ever a negative demand was encountered, the demand series was 
searched backwards for the first positive demand equal to or 
larger than the negative demand. If such a demand was found, 
the positive demand quantity was reduced by the absolute 
value of the negative quantity and the negative quantity was 
set equal to zero. If no offsetting positive demand was 
found, the negative demand quantity was still set equal to 
zero. The resulting edited demand time series was used as 
the base for all further analysis. Two additional demand 
series were created from the edited series: 

1. A demand series of only the non-zero demands, 

2. A demand series of the log of the non-zero demands. 

The sample mean and sample standard deviation for each of the 
three demand series were computed^o and the group, represent- 
ing one line item, was accepted for futher processing if: 



‘ ostationarity was assumed for the demand time series 
under study. If a trend were in fact present, the resampling 
technique effectively would eliminate it by shuffling the 
observed demands into random order. 
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1. There were at least two non-zero demands, and 

2. The standard deviation of all three series was non- 
zero. 

Once accepted, the sample mean and/or sample standard 
deviation of the appropriate series was used to estimate the 
value of Xf such that for a given probability, demand will 
not exceed that value of x. The values of x were computed 
for each model tested at several probabilities from the 
inverse cumulative probability functions. Table 2 displays 
the sample mean and/or sample standard deviation required for 
each model investigated. 

B. THE RESAMPLING PROCEDURE 

The resampling procedure used was chosen over the more 
traditional methods discussed in Chapter II because it pro- 
vides a method of comparing a theotetical distribution to the 
sample data at specified probabilities and is applicable in 
cases such as this where there is a paucity of data available 
for analysis. This procedure is similar to the ’’bootstrap" 
procedure used by Efron [Ref. 8]. The idea behind this pro- 
cedure is to randomly sample with replacement the series of 
available data to create additional pseudo-samples that pos- 
sess the same statistical properties as the original sample. 
The desired statistical property can then be estimated from 
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Table 2: Required Data for Inverse Computations 



flODBL HIM 

1. Exponential A 

2. Normal A 

3. Poisson A 

4. Negative Binomial A 

5. Lognormal C 

6 . Logistic A 

7. La Place A 

8. Bernoulli-exponential 3 

9. Bernoulli- lognorma 1 C 

10. Bernoulli-logistic B 

W here: 



STANDARD 

DEVIAT ION 

N/R 

A 

N/R 

A 

C 

A 

A 

N/R 

C 

B 



A is from edited demand series 

B is from non-zero demand series 

C is from log of non-zero demand series 

N/R indicates parameter Not Required for this model 



each pseudo sample. Wirh a sufficient number of repetitions 
the distribution of the estimated property is known to be 
normal from the central limit theorem with mean equal to the 
population mean and standard deviation equal to the popula- 
tion standard deviation divided by the number of repetitions 
Thus by the use of the resampling procedure, i-: is possi 
ble to study theoretical distributions at points in their 
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right hand tail area, which is the region of most interest in 
computing safety levels in inventory models. The resampling 
procedure, as applied in tais study, is a relatively straight 
forward application of the following technique: 

1. The edited series for each sample of n demand observa- 
tions is treated as the population from which the ran- 
dom samples are drawn . 

2. A random sample of size n is drawn from this popula- 
tion, creating a pseudo sample. The pseudo sample is 
not a permutation of the original population since the 
sample values are selected with replacement from the 
original population. This is easily accomplished in a 
computer by generating n uniform random integers in the 
range from 1 to n, and using these numbers as sub- 
scripts to select the sample from the original 
population. 

3. The pseudo sample percentile is computed for each rep- 
etition from (16) and is compared to the theoretical 
percentile. 

/N 

p = (number of demand observations < x) / n (16) 

The various values of x are computed from the inverse 
cumulative probability functions for each model 
evaluated. 
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Under the null hypothesis that the sample demands come 
from a particular probability distribution, the expected 
value of p should equal p, the theoretical percentile: 

E[ p] = p under Hq (17) 

and the distribution of p should be Normal (p,a2/n) . The per- 
centile p is estimated as above and is estimated from 
standard binomial results as: 

VarCB] = S (1 - p) / n. (18) 



C. MEASURES OF EFFECTIVENESS 

When the above precedure is executed at several theoreti- 
cal percentiles, the fit of the sample data may be compared 
at several locations in the theoretical distribution and the 
overall fit evaluated. Several common measures of effective- 
ness are available to provide a criteria on which to judge 
the success or failure of a model. Some of those available 
are : 

1. Maximum absolute deviation, 

2. Mean squared error, 

3. Algebraic sum of errors. 

Method three was eliminated since errors may offset one 
another resulting in a seemingly good fit, but in reality a 
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very poor fit. Method one is potentially a good way to eval- 
uate error, but is generally not as sensitive as method two. 
Therefore, the mean squared error of the pseudo sample per- 
centile estimates was chosen as the measure of effectiveness 



D. COMPUTER VALIDATION 

The computer programs for this analysis were written in 
FORTRAN IV for execution on the IBM 3033 attached processor 
(System 370) computer installed at the Naval Postgraduate 
School. The programs are batch oriented since the volume of 
data available prevented storage directly on the user’s pri- 
vate disk space. The data tapes received from FMSO were 
loaded onto the system mass storage device and when required, 
all or portions of the data were transferred to the Virtual 
Machine (VM) facility for processing. Two programs were 
written, the first to copy the data to make it accessible to 
the VM system and the second to conducr the actual analysis. 

A complete listing of the FORTRAN source code is available in 
Appendix F. Each program or subprogram contains a documenta- 
tion block at the top that provides a complete description of 



of 'Choice and was accumulated as: 



n 



(19) 




i=l 



34 



the program's purpose, a definition of variables used and a 
list of user-written subroutines and functions required. 

The random number generator used is a part of the 
LLHANDOM II series [Ref. 9] developed at the Naval Postgradu- 
ate School. The programs and subprograms were thoroughly 
tested both independently and as a unit using known data. 

The results were verified by hand held calculator and in all 
cases agreed with the program output. 
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V. 



RESULTS 



The results of the distributional analysis conducted on 
the ten distribution models listed in Table 2 are tabulated 
in Appendix B for 1R consumable items. Appendix B is divided 
into the three demand classes listed in Table 1, and then 
each class is divided by distribution and percentile within 
each distribution. The results for the 2R repairable items 
are presented in a similar format in Appendix C. 

The consistency of the resampling method with varying 
data was evaluated by testing independent sections of the 
consumable and repairable data sets and comparing the 
results. In both cases, independent subsamples were created 
by sampling every third item but changing the initial item 
sampled from the first to the third item. The analysis, not 
presented here, showed that the resampling method produced 
consistent results over all of the subsamples tested. The 
total mean square errcr^^ of the distributions in each sub- 
sample varied by no more than 20^ from the total mean square 
error values obtained from the entire sample population as 
listed in Appendix 3 and Appendix C. The relative standing 



i ‘Total mean square error is the sum of the mean square 
errors computed at the 75th, 80th, 85th, 90th and 95th 
per centiles . 
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q 




of each distribution was not altered in any of the 
subsamples. 

The effect of the number of pseudo sample repetitions was 
also studied. The results proved similar while varying the 
number of repetitions from twenty to fifty. The standard 
deviation of the percentile estimates decreased as the number 
of repetitions increased. The total mean square error for 
each trial varied randomly but remained within 1015 of the 
overall values for each distribution with more repetitions 
producing consistently smaller variations. For the data 
analysis production runs, forty was chosen as the number of 
resampling repetitions as a trade off between computer run 
time and minimizing variance. 

A. LOW-DEMAND ITEMS 

A summary of the total mean squared error figures for the 
models providing the best fit or the smallest total mean 
squared error is provided in Table 3 for the low-demand 
items. The best models for the low demand items all had very 
small total mean squared errors. This can be attributed to 
the very low probability of a demand ever being greater than 
one and to the integer sampling plan for accumulating the 
Bernoulli trials. For the integer results of the Poisson and 
negative binomial d istribut ions, the estimated percentile 
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TABLE 3; 



LOW-DEMAND ITEMS, 

Total Mean Squared Error (10-*) 

Negative 

Po l sson binom ial Conj£Ound 
Consumable 1R 5.64 7.86 10.77 

Repairable 2R 5. 86 7.57 12.91 

will closely match the theoretical percentile unless the num- 
ber of Bernoulli successes is less than the expected number 
of successes for a given theoretical percentile. This pro- 
duces a conservative result that will tend to overstate the 
quantity of stock required by one unit to provide a stated 
level of protection, but in the case of the small quantities 
involved, is not an undesirable preperty. The compound dis- 
tributions, as a group, gave the same total mean squared 
error for much the same reasons. The probability of a non- 
zero demand is very small for the low-demand items, generally 
less than 0.05, thus these models tend to predict a quantity 
of zero for all percentiles. The resampling procedure and 
the sampling plan again provide a conservative result. 

B. MEDIDM-DSMAND ITEMS 

A summary of the total mean squared error figures for the 
models providing the best fit or the smallest total mean 
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squared error is provided in Table 4 for the medium- demand 
items. Here again the predicted demand for the various per- 
centiles is small. The resampling technique still gives the 



TABLE 4: MEDIUM-DEMAND ITEMS, 

Total Mean Squared Error (10“*) 

Negative Bernoulli 

bi n omia l Po isson log no rmal 



Consumable 1R 
Repairable 2R 



28.55 48.28 80.60 

29.14 44.84 112.44 



edge to the integer distributions. The compound distribu- 
tions have begun to show some differences among themselves, 
but as a group are clearly better than any of the continuous 
simple distributions. The magnitude of the total mean 
squared error is significantly larger than for rhe low-demand 
class primarily due to the larger allowable tonal demand 
range defining the medium-demand class. 

C. HIGH-DEMAND ITEMS 

The summary total mean squared error data for the high- 
demand items are listed in Table 5. with the unlimited range 
of mean demands in this class, the higher total mean squared 
errors are expected. The compound models have come into 
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consumable 1R 
Repairable 2R 

their own in this class with the Bernoulli- lognormal yielding 
a total mean squared error less than one half that of the 
normal distribution. The second best distribution, the 
Bernoulli-exponential, gives a total mean squared error only 
slightly higher than the Be rnoulli- lognormal. 



TABLE 5: HIGH-DEMAND ITEMS, 

Total Mean Squared Error (10-'*) 

Bernoulli- Bernoulli- 

l oqno r mal expo nentia l Ex ponenti al 



168. 25 
183. 77 



188.90 

221.80 



298.21 

289.01 
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VI. C ONCLUS ION S AND RE COMM ENDATION S 



A. SUMMARY 

This distributional analysis, conducted using a com- 
pletely different experimental procedure than prior analyses 
provides satisfyingly similar results for the low and medium 
demand categories in showing that the Poisson distribution 
gives the best fit for the low-demand category of items and 
that the negative binomial distribution gives the best fit 
for the medium-demand category of items. The class of com- 
pound distributions, as a whole, gave good fits in the low- 
and medium-demand categories, though the total mean squared 
error was two to four times as large as that of the Poisson 
or negative binomial distribution. 

The analysis indicated that there are several distribu- 
tions which provide better results than the normal distribu- 
tion for high-demand items. The Bernoulli-lognormal 
distribution consistently provided total mean squared errors 
approximately one half that of the normal distribution, indi 
eating a superior fit in the right-hand tail area of the dis 
tribution. Other distributions giving good fits were the 
Bernoulli-exponential and the standard exponential, both 
yielding total mean squared errors less than that of the nor 
mal distribution. 
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B. RECOMMENDATIONS 



The results of this study indicate the Navy should con- 
tinue to use the Poisson and negative binomial distributions 
as the models for the low- and medium-demand classes respec- 
tively. Specifically, the Aviation Supply Office should re- 
evaluate its position on the use of the normal distribution 
for all demand classes and revert to the Poisson and negative 
binomial distributions for low- and medium-demand classes as 
before. The normal distribution tends to inflate the stock 
required for low- and medium-demand items for a specified 
level of protection resulting in increased safety levels and 
excessive dollar investment. 

The Navy should consider replacing the normal distribu- 
tion model used for high-demand items with the compound 
Bernoulli-lognormal distribution. If the Bernoulli-lognormal 
distribution is considered computationally too difficult for 
a large inventory system, the Bernoulli-exponential could be 
used in its place with little loss of effectiveness. 

C. ADDITIONAL RESEARCH 

Follow on research in the following areas may improve and 
expand upon the results presented above: 

1. Optimize the divisions of demand categories (i.e. low, 
medium, high) with the possible use of the probability 
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of a non-zero demand as an element in the 



determination, 

2. Implement the Bernoulli-lognormal model in FMSO's 5A 
simulator to establish the effect of the change on the 
entire supply system, 

3. Investigate the implications of the assumption that the 
demand process is stationary, 

4. Evaluate other compound distributions, such as the 
Bernoulli -log-logisti c, which is particularly appealing 
analytically. 
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AP PSNDIX A 
GLOSS ASY OF TBR MS 

Progra* Related Itea: An item of stock whose demand can be 

predicted from the value of a specific Navy program, i.e., 
flying hours, steaming hours, etc. 

Reorder Point: The on hand stock quantity that when reached, 

triggers an order for replenishment of stock material. The 
reorder Point is the expected demand during procurement lead- 
time plus the safety level. 

Risk: Probability of a stockout during leadtime. 

Safety Level: The quantity of material which is required to 

be on hand to permit continued operation in the event of 
minor interuptions cn normal replenishment or unpredictable 
fluctuations in demand. The safety level determined is 
structured so as to minimize time-weighted, essentiality- 
weighted requisitions short. 

Stockout; A condition that exists when the on hand inventory 
is insufficient to fill the current demand requirements. 
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Systen Perforaance: A subjective measure maximized when the 

time- weighted, essent iality -weighted requisitions short is 
minimized. Time weighting is the consideration of the aver- 
age number of days delay in the availability of material, 
essentiality- weight ing is the consideration of the relative 
essentiality of each item. Requisitions short are requisi- 
tions for which material is not available. 
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A£ PENDIX_B 

CONSaWA BLE ITEM TABOLATED RE SU LTS 
40 REPETITIONS FOR EACH OF 61 PSEODO SAMPLES 



LOW 


DEMAND RANGE; 0.0 


< D < 1 


.0 PER YEAR 


DISTRIBUTION; 


POISSON 














THEORETICAL PERCENTILES (P) 




ESTIMATES; 


Oi.25 


o| 

0O| 

•1 

Ol 


0 2.85 


0.90 


0^95 


MEAN 


0. 7500 00 


0.81 25 00 


0.854006 


0.913924 


0.956069 


STD DEV 


0. 008766 


0.007902 


0.007148 


0.005678 


0.004149 


MSE 


o 

• 

o 


0.00 0156 


0.000018 


0.000290 


0.000100 


TOTAL MSE; 


5. 64E- 04 










DISTRIBUTION; 


NEGATI VE 


I BINOMIAL 












THEORETICAL PERCENTILES (P) 




ESTIMATES; 




0^80 


Q-.85 


0. 90 


0.95 


MEAN 


0. 750000 


0.31 2500 


0.854006 


0.913924 


0. 948994 


STD DEV 


0. 008766 


0 .007902 


0.007148 


0.005678 


0. 004454 


MSE 


o 

• 

o 


0.000156 


0.000018 


0. 000290 


0.000322 



TOTAL MSE: 7.86E-04 



46 



40 REPETITIONS FOR EACH OF 61 PSEUDO SAMPLES 



LOW 


DEMAND 


RANGE: 0.0 


< D < 1. 


0 PER YEAR 




DISTRIBUTION: 


NORMAL 














THEORETICAL PERCENTILES (P) 




E^IMAOS: 


0.75 


0^80 


Oil 5 


Oiio 


0.95 


MEAN 


0.953221 


0.95 3221 


0.953221 


0.953221 


0.953221 


STD DEV 


0. 004275 


0.004275 


0.004275 


0.004275 


0.004275 


MSE 


0. 0423 19 


0.024488 


0.01 1662 


0.003835 


0.001006 


TOTAL MSE: 


833. 09E- 


04 









DISTRIBUTION: EXPO NENTI AL 







THEORETICAL PERCENTILES (P) 




E^IMMES: 


0i25 


Oj.80 


0^85 


OiiO 


0i95 


MEAN 


0. 953221 


0.953221 


0.953221 


0.953221 


0.953221 


STD DEV 


0. 004275 


0.004275 


0.004275 


0.004275 


0.004275 


MSE 


0. 0423 19 


0.024488 


0.01 1662 


0. 003835 


0.001006 


TOTAL MSE: 


833. 09E- 


04 









DISTRIBUTION; LOGNORMAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 




0^0 


2^85 


0^22 


2^.95 


MEAN 


0. 9784 04 


0.995246 


0.995246 


1.000000 


1.000000 


STD DEV 


0. 002943 


0.00 1392 


0.001392 


O 

• 

o 


o 

• 

o 


MSE 


0. 052629 


0.03 8292 


0.021267 


0.009998 


0. 002500 



TOTAL MSE: 1246.85E-04 
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40 REPETITIONS FOR EACH OF 61 PSEUDO SAMPLES 
LOW DEMAND RANGE; 0.0 < D < 1.0 PER YEAR 



DISTRIBUTION; LAPLACE 



ESTIM ATES ; 
MEAN 
STD DEV 
MSE 

TOTAL MSE; 



THEORETICAL PERCENTILES (P) 





•o 

1*00 

lo 


0i85 


0. 9 0 


0^95 


0.953221 


0.953221 


0.953221 


0.953221 


0.953221 


0. 004275 


0.004275 


0.004275 


0.004275 


0.004275 


0. 042319 


0.024488 


0.01 1662 


0.003835 


0.001006 



833. 09E-04 



DISTRIBUTION; LOGISTIC 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


Q.2.11 




2i85 


0i20 


0.95 


MEAN 


0.953221 


0.953221 


0.953221 


0.953221 


0.953221 


STD DEV 


0. 004275 


0.004275 


0.004275 


0.004275 


0.004275 


MSE 


0. 042319 


0.024488 


0.01 1662 


0.003835 


0.001006 



TOTAL MSE; 833.09E-04 



DISTRIBUTION 

ESTIMATES; 
MEAN 
STD DEV 
MSE 

TOTAL MSE; 



BER NOU LLI :iEX PONE NTIii 

THEORETICAL PERCENTILES (P) 



Oi.25 

0. 7500 00 
0. 008766 
0 . 0 

10.77S-04 



2i80 

0.81 2500 
0.007902 
0.000156 



2i65 

0.854006 

0.007148 

0.000018 



0.90 

0.913924 
0. 005678 
0. 000290 



0 .95 

0. 944939 
0.004618 
0.000613 
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40 REPETITIONS FOR EACH OF 61 PSEUDO SAMPLES 



LOH DEMAND RANGE; 0.0 < D < 1.0 PER YEAR 



DISTRIBUTION; BSRMOO LLI -LO GNORM AL 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


0i25 


0^90 




OiiO 


0^95 


MEAN 


0. 750000 


0.31 25 00 


0.854006 


0.913924 


0.944939 


STD DEV 


0. 008766 


0.007902 


0.007148 


0.005678 


0.004618 


MSE 


O 

• 

o 


0.000156 


0.000018 


0.000290 


0.000613 



TOTAL MSE: 10.77E-04 

DISTRIBUTION; BERNOULLI-LOGISTIC 



ESTIM ATE S; 
MEAN 
STD DEV 
MSE 

TOTAL MSE; 



0.75 

0. 750000 
0. 008766 
0 . 0 



THEORETICAL PERCENTILES (P) 



0^80 
0.81 2500 
0.00 79 02 
0.000156 



0^85 

0.354006 

0.007148 

0.000018 



OiiO 

0. 913924 
0.005678 
0.000290 



Oils 

0.944939 

0.004618 

0.000613 



10.77E-04 
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40 REPETITIONS FOR EACH OF 753 PSEODO SAMPLES 



MEDIUM 


DEMAND RANGE: 1.0 


< D < 20.0 


PER YEAR 




DISTRIBUTION: 


POISSON 














THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oi.25 


0 .80 


0^.85 


0-.20 


0.95 


MEAN 


0. 7481 31 


0.807719 


0. 843932 


0.893923 


0. 925303 


STD DEV 


0. 002501 


0.00 2271 


0. 002091 


0.001774 


0.001515 


MSE 


0. 0001 38 


0.00 0398 


0. 000603 


0.001423 


0.002266 


TOTAL MSE: 


48. 28E-04 










DISTRIBUTION: 


NEGATIVE 


_ BINOMIAL 












THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0^75 


0.1 80 


8^85 


OiiO 




MEAN 


0. 746099 


0.80 5430 


0.843400 


0.900210 


0.943628 


STD DEV 


0. 002508 


0.002281 


0.002094 


0.001727 


0.001329 


MSE 


0. 000332 


0.0005 22 


0.000589 


0.000807 


0.000606 


TOTAL MSE: 


28.55E-04 










DISTRIBUTION: 


NORMAL 














THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oils 


0^80 




Oi-lO 


0.95 


MEAN 


0. 888030 


0.896699 


0.907882 


0.920152 


0.935573 


STD DEV 


0. 001817 


0.00 1754 


0. 00 1666 


0.001562 


0.001415 


MSE 


0. 024771 


0.01 4168 


0.007327 


0.003509 


0.002400 



TOTAL MSE: 521.74E-04 
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40 REPETITIONS FOR EACH OF 753 PSEDDO SAMPLES 



MEDIUM DEMAND RANGE; 1.0 < D < 20.0 PER YEAR 



DISTRIBUTION: EXP ONE NTIAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


Oils 


0 .80 


Qi85 


0i90 


0.95 


MEAN 


0. 852981 


0.86 1003 


0.871591 


0.884311 


0.898352 


STD DEV 


0. 002040 


0.00 1993 


0. 001928 


0.001843 


0.001741 


MSS 


0, 020357 


0.01 1779 


0.006819 


0.004953 


0.006210 



TOTAL MSE; 501.19E-04 



DISTRIBUTION; LOGNORMAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


Oils 


2i80 


0i85 


0.90 


0.95 


MEAN 


0.953018 


0.96 2552 


0.970638 


0.980623 


0. 990731 


STD DEV 


0. 001219 


0.00 1094 


0.000973 


0.000794 


0.000552 


MSE 


0, 043846 


0.02 8367 


0.01 5977 


0.007411 


0.002040 



TOTAL MSE; 976.41E-04 



DISTRIBUTION: LAPL ACE 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


0^75 


Qi80 


2i85 


2i20 


0.95 


MEAN 


0. 884254 


0.893767 


0.904733 


0.9181 10 


0.936078 


STD DEV 


0, 001843 


0.00 1775 


0.00 1692 


0.001580 


0.001409 


MSE 


0. 024079 


0.01 3865 


0.007191 


0.003573 


0.002363 


TOTAL MSE; 


510. 72E- 


04 
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40 REPETITIONS FOR EACH OF 753 PSEUDO SAMPLES 
MEDIUM DEMAND RANGE: 1.0 < D < 20.0 PER YEAR 



DISTRIBUTION: LO GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


0.1 80 


OiM 


0^.90 


0^95 


MEAN 


0. 875036 


0.886687 


0.896908 


0.913037 


0.934834 


STD DEV 


0, 001905 


0.00 1826 


0.001752 


0.001624 


0.001422 


MSE 


0. 022871 


0.01 3193 


0.006955 


0.003690 


0.002447 


TOTAL MSE: 


491 . 56E- 


04 









DISTRIBUTION: BE HNOO LLI-EX PONENTIAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 




Qi80 


0.85 


0 . 2 . 10 . 


0.95 


MEAN 


0. 738072 


0.795951 


0.840430 


0.896563 


0.950635 


SID DEV 


0. 002533 


0.002322 


0. 0021 10 


0.001755 


0.001248 


MSE 


0. 002132 


0.00 2338 


0.002267 


0.002224 


0.001327 


TOTAL MSE: 


102.88E- 


04 








DISTRIBUTION: 


3ERN0U LLI-LOGNORMAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 




Qi80 


2 1.8 5 






MEAN 


0. 749344 


0.809159 


0. 849391 


0.903886 


0.948263 


STD DEV 


0. 002497 


0.002264 


0.002061 


0.001698 


0.001276 


MSE 


0. 001456 


0.00 1751 


0.001730 


0.001866 


0.001253 



TOTAL MSE: 30.60E-04 
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40 REPETITIONS FOR EACH OF 753 PSEODO SAMPLES 



MEDIUM DEMAND RANGE: 1.0 < D < 20.0 PER YEAR 



DISTRIBUTION: BERNOU LLI-L O GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 




0^80 






0.95 


MEAN 


0. 7600 06 


0.82 1542 


0.864176 


0.916219 


0.955416 


STD DE7 


0. 002461 


0.002206 


0.001974 


0.001596 


0.001189 


MSB 


0. 002256 


0.00 24 99 


0. 002304 


0.002030 


0.001241 


TOTAL MSE: 


103. 30E- 


04 
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40 REPETITIONS FOR EACH OF 429 PSEUDO SAMPLES 



HIGH 


DEMAND RANGE: 20.0 


< D PER 


YEAR 




DISTRIBUTION 


: POISSON 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0^75 0^0 




0.90 


0^95 


MEAN 


0.689018 0.730300 


0.753320 


0.780038 


0. 800697 


STD DEV 


0. 003534 0.00 3388 


0. 003291 


0.003162 


0.003050 


MSE 


0.024843 0.027146 


0.033935 


0. 0421 54 


0. 051388 


TOTAL MSE: 


1794. 66E-04 








DISTRIBUTION 


: NEGATIVE BINOMIAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0^7 5 0^0 




0-.90 




MEAN 


0.724173 0.747130 


0.7913 13 


0.848041 


0.896526 


STD DEV 


0. 0034 12 0.00 3318 


0. 0031 02 


0.002740 


0. 002325 


MSE 


0.003047 0.024647 


0.026816 


0,029773 


0.031637 


TOTAL MSE: 


1 159. 20E-04 








DISTRIBUTION 


: NORMAL 










THEORETICAL PERCENTILES (P) 




estimates: 


0. 75 0^0 


2i.85 


0^.90 


0.95 


MEAN 


0. 850264 0.87 1100 


0.890154 


0.910670 


0.934248 


STD DEV 


0.002724 0.002558 


0. 002387 


0.002177 


0.001892 


MSE 


0. 016622 0.01 0312 


0.005815 


0.003155 


0.002227 



TOTAL MSE: 38 1. 32 E- 04 
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40 HEPETITIONS FO H BACH OF 429 PSEODO SAMPLES 



HIGH 


DEMAND 


RANGE: 20.0 


< D PER 


YEAR 




DIS TRIBUTION: 


exponenti^ 












THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oill 


0^8 0 


0^.85 


OiiO 


0.95 


MEAN 


0. 784610 


0.81 4275 


0.846685 


0.881207 


0.918702 


STD DEV 


0. 003138 


0.00 2969 


0. 002750 


0.002470 


0.002086 


MSE 


0. 009534 


0.006885 


0.005219 


0.004384 


0.003800 


TOTAL MSE: 


29 8. 2 IE- 


04 









DISTRIBUTION: LOGNORMAL 







THEORETICAL PERCENTILES (P) 




ESTIMATES: 




0^80 


8.85 




0.95 


MEAN 


0. 855867 


0.886448 


0.918614 


0.949144 


0.978288 


STD DEV 


0- 002681 


0.00 2422 


0. 002087 


0.001677 


0.001113 


MSE 


0. 019028 


0.01 3031 


0.0081 13 


0.004389 


0.001639 



TOTAL MSE; 462.00E-04 



DISTRIBUTION: LAPL ACE 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^25 


QiiO 


8.85 


8i20 


0.95 


MEAN 


0. 843416 


0.86 3600 


0.885612 


0. 907473 


0.934956 


STD DEV 


0. 002774 


0.00 2620 


0.002430 


0.002212 


0.001883 


MSE 


0. 015775 


0.00 9709 


0.005680 


0.003248 


0. 002167 


TOTAL MSE: 


365.79E- 


04 
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40 REPETITIONS FOR EACH OF 429 PSEUDO SAMPLES 



HIGH DEMAND RANGE: 20.0 < D PER YEAH 



DISTRIBUTION: LOGISTIC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 




0^90 


2i85 


0i20 




MEAN 


0. 822132 


0.84 6791 


0.871709 


0.9001 10 


0.933755 


STD DEV 


0. 002919 


0.002750 


0.002553 


0.002289 


0. 001899 


MSE 


0. 013721 


0.008841 


0.005633 


0.003523 


0. 002259 



TOTAL MSE: 339.76E-04 



DISTRIBUTION: BER NOULLI-E X PONENTIAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


0^0 


2*.85 


Oi.10 


0^.95 


MEAN 


0. 759754 


0.81 6600 


0. 867492 


0.914225 


0.955101 


STD DEV 


0. 003261 


0.002954 


0.002588 


0. 002138 


0. 001531 


MSB 


0. 005828 


0.00 5034 


0.003977 


0.0026 18 


0. 00 1432 



TOTAL MSE: 188.90E-04 



DISTRIBUTION: BS HNOD LLI -L O GNORM AL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


0 .80 


Os.85 






MEAN 


0. 745078 


0.797614 


0. 847684 


0.901963 


0.953670 


STD DEV 


0. 003327 


0.003067 


0.002743 


0. 002270 


0. 001605 


MSE 


0. 004768 


0.004403 


0.003632 


0.0026 13 


0. 00 1409 



TOTAL MSE: 168.25E-04 



56 



40 REPETITIONS FOR EACH OF 429 FSEDDC SAMPLES 



HIGH DEMAND RANGE: 20.0 < D PER YEAR 



DISTRIBUTION: BERNOULLI -LOGISTIC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0.75 


0.80 


Qi.85 


0i90 


0. 95 


MEAN 


0.828539 


0.86 4722 


0.892571 


0.921326 


0 .947036 


STD DEV 


0. 002877 


0.00 2611 


0.002364 


0. 002055 


0.001710 


MSE 


0. 012517 


0.00 8862 


0.005485 


0.003020 


0.001538 



TOTAL MSE: 314.23E-04 
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AP PENDI X C 

REP AIRABLE IT EM.! ABOL AT ED_R S SUITS 
40 REPETITIONS FOR EACH OF 137 PSEUDO SAMPLES 
LOW DEMAND RANGE; 0.0 < D < 1,0 PER YEAR 



DISTRIBUTION; POISON 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


0^75 


0.80 




OiiO 


0.95 


MEAN 


0. 7500 00 


0.81 2496 


0.353613 


0.912594 


0. 955936 


STD DEV 


0. 005849 


0.005273 


0.004775 


0.003815 


0.002772 


MSE 


0. 0 


0.000156 


0.000019 


0.000311 


0,000099 



TOTAL MSE; 5.86E-04 



DISTRIBUTION: NEGATIVE BINOMIAL 



ESTIMATES; 
MEAN 
STD DEV 
MSE 

TOTAL MSE: 



THEORETICAL PERCENTILES (P) 
0.80 



0^11 
0. 7500 00 
0. 005849 
0.0 

7.57E-04 



0.81 2496 
0.005273 
0.000156 



0.853613 

0.004775 

0.000019 



OiiO 

0.912594 

0.003815 

0.000311 



0.95 

0.950753 

0.002923 

0.000270 
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40 REPETITIONS FOR EACH OF 137 PSEUDO SAMPLES 



LOW 


DEMAND 


RANGE; 0.0 


< D < 1 


. 0 PER YEAR 




DISTRIBUTION; 


NORMAL 














THEORETICAL PERCENTILES (P) 




ESTIMATES; 


0.75 


0.80 








MEAN 


0. 949164 


0.94 9164 


0.949164 


0.949 164 


0.949164 


STD DEV 


0. 002967 


0.002967 


0.002967 


0.002967 


0.002967 


MSE 


0. 040918 


0.023465 


0.01 1016 


0.003568 


0. 001118 


TOTAL MSE; 


800. 84E- 


04 









DISTRIBUTION: E XPO NENTIAL 









THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0.75 




0.80 




OiiO 


0.95 


MEAN 


0. 949164 


0 


.94 9164 


0. 949 1 64 


0.949.164 


0.949164 


STD DEV 


0. 002967 


0 


.00 2967 


0.002967 


0.002967 


0.002967 


MSB 


0. 040918 


0 


.02 3465 


0.01 1016 


0.003568 


0.001118 


TOTAL MSE: 


800. 84E- 


04 











DISTRIBUTION; LpGSOR^L 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 






0^85 




0. 95 


MEAN 


0. 978874 


0.99 1028 


0. 99 1028 


1.000000 


1. 000000 


STD DEV 


0. 001943 


0.00 1274 


0.00 1274 


o 

• 

o 


O 

• 

o 


MSE 


0. 052955 


0.036815 


0. 020201 


0.009996 


0.002500 



TOTAL MSE: 1224.66E-04 



59 



40 REPETITIONS FOR EACH OF 137 PSEODO SAMPLES 



LOW DEMAND RANGE; 0.0 < D < 1.0 PER YEAR 



DISTRIBUTION; LAPLACE 



ESTIM ATES : 
MEAN 
STD DEV 
MSB 

TOTAL MSE; 



THEORETICAL PERCENTILES (P) 



0.75 


0.80 


0.85 


0.90 


0.95 


0.949164 


0.94 9164 


0.949164 


0.949164 


0.949164 


0. 002967 


0.002967 


0.002967 


0.002967 


0.002967 


0. 0409 18 


0.023465 


0. 01 10 16 


0. 003568 


0.001118 



800. 84E-04 



DISTRIBUTION; LO GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES; 


Oils 


0i80 


OiM 


OiiO 


0.95 


MEAN 


0. 9491 64 


0.94 9164 


0.949164 


0.949164 


0.949164 


STD DEV 


0. 002967 


0.002967 


0.002967 


0.002967 


0.002967 


MSE 


0. 040918 


0.02 3465 


0.01 1016 


0.003563 


0.001113 


TOTAL MSE: 


800. 84E- 


04 









DISTRIBUTION; BERNOULLI-EX PONENTIAL 









THEORETICAL PERCENTILES ( P) 




ESTIMATES: 


Oils 




OiSO 


Oil 5 


0.9 0 


0i95 


MEAN 


0. 750000 


0 


.81 2496 


0.853613 


0.912594 


0.942635 


STD DEV 


0. 005849 


0 


.005273 


0.004775 


0.003815 


0.003141 


MSE 


O 

• 

o 


0 


.00 0156 


0.000019 


0.000311 


0.000804 


TOTAL MSE; 


12. 91E- 


04 
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40 REPETITIONS FOR BACH OF 137 PSEDDO SAMPLES 



LOW 


DEMAND 


RANGE: 0.0 


< D < 1 


.0 PER YEAR 




DISTRIBUTION: BERNOULLI-LOGNORMAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0-25 


oi 

OI 


0^85 


0-.20 


0.95 


MEAN 


0. 750000 


0.81 2496 


0.853613 


0.912594 


0.942635 


STD DEV 


0. 005849 


0.005273 


0.004775 


0.003815 


0.003141 


MSE 


O 

• 

o 


0.00 0156 


0.000019 


0.000311 


0.000804 


TOTAL MSE: 


12.91E- 


04 









DISTRIBUTION: BERNOUL LI-L OGISTIC 



ESTIM ATES : 
MEAN 
STD DEV 
MSB 

TOTAL MSE; 



0.75 

0. 7500 00 
0. 005849 
0 . 0 



THEORETICAL PERCENTILES (P) 



0^80 
0.81 2496 
0.005273 
0.000156 



2i85 

0.853613 

0.004775 

0.000019 



0^20 

0.912594 

0.003815 

0.000311 



0^85 

0.942635 
0.003141 
0. 000804 



12. 91E-04 
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aO REPETITIONS FOR BACH OF 1487 PSEUDO SAMPLES 



MEDIUM 


DEMAND 


RANGE: 1.0 


< D < 20.0 


PER YEAR 




DISTRIBUTION: 


POISSO 


N 












THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0.75 


0^80 




Oi.90 


0.95 


MEAN 


0. 747034 


0.806323 


0.843928 


0.894732 


0.929457 


STD DEV 


0. 001782 


0.00 1620 


0.001488 


0. 001258 


0.001050 


MSE 


0.000238 


0.000507 


0.000582 


0.001321 


0.001836 



TOTAL HSE: 44.84E-04 



DISTRIBUTION: NEGATIV E BI NOMIAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


0^80 






0^85 


MEAN 


0. 744688 


0.304450 


0.843869 


0.900645 


0.944953 


STD DEV 


0. 001788 


0.00 1626 


0.001488 


0.001227 


0.000935 


MSE 


0. 000444 


0.00 0598 


0.000533 


0.000788 


0.000551 



TOTAL MSE; 29.14E-04 



DISTRIBUTION: NORM^ 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


lO 

1*00 

lo 


2i85 


0-.20 




MEAN 


0.861474 


0.87 5334 


0.890316 


0.905783 


0.927155 


STD DEV 


0. 001416 


0.00 1354 


0.00 1 281 


0. 001198 


0.001066 


MSE 


0. 019444 


0.01 1363 


0.006084 


0.003600 


0.002944 


TOTAL MSE: 


434. 35E- 


04 
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40 REPETITIONS FOR EACH OF 1487 PSEUDO SAMPLES 



MEDIUM DEMAND 


RANGE: 1.0 


< D < 20.0 


PER YEAR 




DISTRIBUTION: EXPONENTIAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 0^75 


2^80 






0^95 


MEAN 0.823217 


0.834022 


0.845365 


0. 86427 1 


0.886389 


STD DEV 0.001564 


0.00 1526 


0.001482 


0. 001404 


0.001301 


MSB 0.016435 


0.01 0140 


0.007542 


0.006691 


0.007993 


TOTAL MSE: 488.02E- 


04 









DISTRIBUTION: LOG NORMAL 







THEORETICAL PERCENTILES (P) 




ESTIMATES: 




0.80 


2i85 




0.95 


MEAN 


0. 946618 


0.955351 


0.964288 


0. 975138 


0.987135 


STD DEV 


0.000922 


0.000847 


0.000761 


0.000638 


0.000462 


MSB 


0. 041761 


0.026454 


0. 01 4803 


0.006692 


0.001890 


TOTAL MSE: 


916. OOE- 


04 









DISTRIBUTION: LAPLACE 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^25 


0 .80 


22.85 


22.90 


0.95 


MEAN 


0. 856232 


0.87 0161 


0.885855 


0. 902481 


0.928262 


STD DEV 


0.001439 


0.00 1378 


0. 00 1304 


0.001216 


0.001058 


MSE 


0. 018823 


0.01 1026 


0.006113 


0.003736 


0. 002851 



TOTAL MSE: 425.48E-04 
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40 REPETITIONS FOR EACH OF 1487 PSEUDO SAMPLES 



MEDIUM 


1 DEMAND 


RANGE: 1.0 


< D < 20.0 


PER TEAR 




DISTRIBUTION: 


LOGISTIC 












THEORETICAL PERCENTILES (P) 




ESTIMATES; 


Oils 


0^80 


0^85 


0i80 


Oils 


MEAN 


0.845409 


0.858922 


0.875792 


0.896303 


0.925823 


STD DEV 


0.001482 


0.00 1427 


0.001352 


0.001250 


0.00 1075 


MSE 


0. 017690 


0.01 0603 


0.006269 


0.004086 


0.003093 


TOTAL MSE; 


417.42E- 


04 








DISTRIBUTION: 


BERNOULLI-EX PONENTIAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oils 


0^80 


Qi85 




OiiS 


MEAN 


0.730065 


0.787917 


0.834508 


0.893026 


0.954267 


STD DEV 


0. 001820 


0.00 1676 


0.00 1524 


0.001267 


0.000857 


MSB 


0. 003437 


0.00 3608 


0.003449 


0.003021 


0. 001465 


TOTAL MSE: 


149.79E- 


04 








DISTRIBUTION: 


BE^^IJ;I:^g FORMAL 







THEORETICAL PERCENTILES (P) 



ESTIMATES: 


Oils 


0.80 


0^85 


OiiO 


0^95 


MEAN 


0. 750179 


0.80 8666 


0.852134 


0. 904194 


0. 949477 


STD DEV 


0. 001775 


0.00 1613 


0.00 1455 


0.001207 


0.000898 


MSE 


0.002226 


0.00 2681 


0.002600 


0.002327 


0.001410 



TOTAL MSE: 112.44E-04 
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40 REPETITIONS FOR EACH OF 1487 PSEODO SAMPLES 



MEDIUM DEMAND RANGE: 1.0 < D < 20.0 PER YEAR 



DISTRIBUTION: BE RNOU LLI-LO GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0.75 


0^8 0 


2 s. a 5 


2s.M 


0.95 


MEAN 


0.762406 


0.82 2694 


0.866844 


0.915851 


0,953988 


STD DEV 


0. 001745 


0.00 1566 


0,00 1393 


0.001138 


0.000859 


MSB 


0. 003183 


0.003523 


0. 0031 07 


0.002468 


0.001379 



TOTAL MSE: 136.61E-04 
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40 REPETITIONS FOR EACH OF 356 PSEODO SAMPLES 



HIGH 


DEMAND RANGE: 20.0 


< D PER 


YEAR 




DISTRIBUTION 


: POISSON 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0.75 0.80 


0i85 


0^90 


Oils 


MEAN 


0.695829 0.735554 


0.762330 


0.791447 


0. 821025 


STD DEV 


0. 003855 0.00 3696 


0.003567 


0.003405 


0. 003212 


MSB 


0.022812 0.026569 


0.032141 


0. 038638 


0.044655 


TOTAL MSE: 


1648. 16E-04 








DISTRIBUTION 


: NEGATIVE BINOMIAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oils 0.80 


OiM 


0i90 


Oils 


MEAN 


0.722861 0.755975 


0.799852 


0. 857560 


0.904930 


STD DEV 


0.003751 0.003599 


0. 003353 


0.002929 


0.002458 


MSE 


0.005285 0.024724 


0.026795 


0.028699 


0.030113 


TOTAL MSE: 


1 156. 16E-04 








DISTRIBUTION 


: NORMAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


Oils Oj.80 


Oils 


OiiO 


Oils 


MEAN 


0.836334 0.859810 


0.883941 


0.905503 


0.932035 


STD' DEV 


0. 0031 00 0.00 2909 


0. 002684 


0.002451 


0. 002109 


MSE 


0.014669 0.009278 


0.005429 


0.003190 


0. 002319 


TOTAL MSE: 


348. 85E- 04 
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40 REPETITIONS FOR EACH OF 356 PSEODO SAMPLES 



HIGH DEMAND RANGE: 20.0 < D PER YEAR 



DISTRIBOTION: EXPONEN^AL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


Oils 


0.80 


Qi85 


lo 

lo 


0.95 


MEAN 


0.780940 


0.81 8187 


0. 856974 


0. 898153 


0.935255 


STD DEV 


0. 003466 


0.003232 


0.002934 


0. 002535 


0.002062 


MSB 


0. 009578 


0.00 7296 


0. 005229 


0.003830 


0. 002968 



TOTAL HSE: 289.01E-04 



DISTRIBUTION: LOGNOR MAL 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0.75 


0^0 


0i85 


OiiO 


0.95 


MEAN 


0. 8191 19 


0.855337 


0.895522 


0.933281 


0.971445 


STD DEV 


0. 003226 


0.002948 


0. 002563 


0.002091 


0.001396 


MSE 


0. 013003 


0.00 8704 


0. 005985 


0. 003420 


0.001454 



TOTAL MSE: 325. 65E-04 



DISTRIBUTION: LAPL ACE 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 




0^80 


0^85 


Oi^O 




MEAN 


0. 827652 


0 .852088 


0.877334 


0. 902360 


0.932870 


STD DEV 


0. 003165 


0.00 2975 


0.002749 


0. 002487 


0.002097 


MSE 


0. 013659 


0.008798 


0. 005346 


0.003302 


0.002267 


TOTAL MSE: 


333. 71E- 


04 
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40 REPETITIONS FOR EACH OF 356 PSEODO SAMPLES 



HIGH DEMAND RANGE; 20.0 < D PER YEAR 



DISTRIBOTION: LO GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


0^75 


0^0 


02.85 


0i£0 


£ 2.25 


MEAN 


0.805812 


0.83 3134 


0.860882 


0.892308 


0. 930748 


STD DEV 


0. 003315 


0.003125 


0.002900 


0.002598 


0.002128 


MSE 


0. 012232 


0.008347 


0.005666 


0. 003803 


0.002397 


TOTAL MSE; 


324. 44E- 


04 








DISTRIBUTION: 


BERNODLLI-EX RONE NT 


'Ikh. 










THEORETICAL PERCENTILES (P) 




ESTIMATES; 


0iI5 


0^90 


Q 2.85 


O 2 . 8 O 


0.95 


MEAN 


0. 770942 


0.829521 


0. 877192 


0.921468 


0.S60362 


STD DEV 


0. 003522 


0.00 3151 


0.002750 


0.002254 


0. 00 1635 


MSE 


0. 0071 08 


0.005916 


0.004548 


0.003092 


0.001516 


TOTAL MSE; 


221.80E- 


04 








DISTRIBOTION; 


BERNOULLI-LOGNORMAL 










THEORETICAL PERCENTILES (P) 




ESTIMATES: 


0^75 


0.80 


O 1.85 


O 2 . 9 O 


0.95 


MEAN 


0. 7461 08 


0.797150 


0. 847197 


0.901039 


0.953617 


STD DEV 


0. 003647 


0.00 3370 


0.003015 


0.002502 


0.001762 


MSB 


0. 005655 


0.00 4764 


0.003680 


0.002809 


0.00 1469 


TOTAL MSE: 


183.77E- 


04 
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UO REPETITIONS FOR EACH OF 356 PSEDDO SAMPLES 



HIGH DEMAND RANGE: 20.0 < D PER YEAR 



DISTRIBUTION: BSRHOOLLI-LO GIST IC 

THEORETICAL PERCENTILES (P) 



ESTIMATES: 


Oill 




2 =.3 5 


0^90 


0^95 


MEAN 


0. 829834 


0.86 1163 


0.886755 


0.914062 


0.94 1856 


STD DEV 


0.003149 


0.00 2898 


0.002656 


0. 002349 


0.001961 


MSE 


0. 013074 


0.00 8963 


0.005495 


0.003143 


0.001722 



TOTAL MSE: 323.98E-04 
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AP PENDIX D 



ESTIMATION CF LOGNORMAL MEA N AND VA RIAN CE FROM U NTR ANSFORMED 

DAT A 



The usual method of estimating the sample mean and vari- 
ance of a random sample assumed to come from a lognormal dis- 
tribution is to take the log transform of the data and then 
compute the standard method of moments estimates. This is 
time consuming on a computer as the logarithm must be taken 
for each observation in the sample. Formulae xo compute the 
mean and variance of the un transformed data from the mean and 
variance of the transformed data can be derived for the log- 
normal distribution. 

If Y is distributed nor mal(Uy ,Oy) , then X is distributed 
lognormal if Y = ln(X). The mean of X [Ref. 10] is expressed 
as: 



Ux = 





( 20 ) 



and the variance as; 





( 21 ) 



Equations (20) and (2 1) are the reverse of that required in 
this case since and are obtained from the transformed 
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data. What is needed then, is to solve for \iy and in 
teras of and a|. Substituting (20) into (21) results in: 




aZ 

which when solved for e 
the solution of op as: 



( 22 ) 

talcing logs and simplifying, gives 




(23) 



Substituting (23) back into the log of (20) , solving for Uy 
and simplifying gives the solution of py as: 



(2U) 




Thus the mean and variance cf the transformed data may be 
estimated from the sample mean and variance of the untrans- 
fcrmed data reducing significantly the computational 
requirements. 
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CONDUCT RESAMPLING PROCEDURE RETURNING THE NUMBER OF 
BERNOULLI SUCCESSES AND THE AVERAGE MEAN SQUARED ERROR OF THE 
MODEL. 
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END 
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EXTRACT CHARACTERS ONE AT A TIME FROM IVAL AND STORE AS 
ALPHANUMERIC CHARACTERS IN ARRAY NUM 
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CREATE OUTPUT ARR AY 



STR$ = B$ 

IF (IPRT .GE. 30) WRITE (ICNOUT, 220) IVAL,STR$ 

220 FORMATC NUMERIC VALUE = ',18,' ALPHANUMERIC VALUE 
RETURN 
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NBRNEQ — COUNTS NUMBER OF SAMPLES WITH NEGATIVE OBS. 
NBROBS -- COUNTS OF OBSERVATIONS IN EACH OF ABOVE 

NEGOBS — LOGICAL VARIABLE, TRUE IF SAMPLE HAS NEG OBS 

0BS(H,1) -- ORIGINAL OBSERVATIONS WITH NEGATIVES REMOVED 

C*;-^ 0BS(M,2) — ALL ZERO OBSERVATIONS REMOVED FROM DATA 
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IF (IPRT ,GE. 5) WRITE (IPRINT, 630) (OBS (K, 3) , K= 1 , CO UNT) 
RETURN 

40 CONTINUE 
RETURN 1 

600 FORMATC 1RAW OBSERVATIONS • , 3 (T20, 20F5 . 0 ,/) ) 



610 FORMATCONEGATIVB REMOVED • , 3 (T20, 20F5. 0,/) ) 

620 FORMAT('0NON ZERO OBS ONLY' , 3 (T20, 20F5. 0,/) ) 

630 FORMATCOLOG NON-ZERO OHL Y • , 3 (T20 , 20F5 . 0 ,/) ) 
800 FORMAT (50F4 .0) 

END 
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DEM(1550), MONTH(NBR), - QUANT(46) 

LOG ICAL>^1 
TIME 

COMMON /FILES/ 
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DO 80 J=FIRST, LAST 
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INVERSE STANDARD NORMAL FUNCTION 



C>3 



W 

cu 

>4 

EH 

s 



X 

ac 

o 



tsl 



M 

CU 

X 

£h 

> 

Q 

EH 

CO 



W 

cu 

>H 

H 

E 

> 

O 

EH 

cn 





a 






0 


0 






25 




-f 








pq 












Eh 






M 


K-} 






0 












cu 












CO 






Eh 








M 










z 


X 


















U 


•» 






Eh 


^0 








0 


£h 


0 






z 


« m 


+ 






2 








U 


03 








M 


E 


:a* 






0 


W 








O 








E 


W T- 








£h 




r— 






H 


H 








b 








E 


H 


W 






U 


E 








EH 


^ O 










rD 






P^4 


0 


04 






Z 


E 


0 






U 




rD 






I-:; 










EH 


X 






E 


% 


EH 






Z 










-S2 










0 


Eh 






0^ 




0 






E 




M 






M 


cu 






< 


^ 0 


S 
















Pm 


u:: 


CU 






EH 


X 






E 


u: 








CJ 


% 


03 










X 






2S 


Eh 






« 




E 






M 


E 


pq ^ 






pq 


X o 


Eh 






W 


E 






n 


X 0 


E 






Eh 




H 0 






CJ 


Ct4 


S 






Z 








25 










CO 


X 








C 


t 








0 


S 








^ • 


Oi 






H 




t 






hj 


cn w 








CX 4 


E 






0 


CO w 


X 






E 


EH 


li^; pq 






Pl4 


H 2 : 


s 






X 


1 






0 


M 25 


w 






0 


CO 








<3 


cn t 








w 










03 « 








hO 


H 


X • 






h3 




II 








II 








E 


II 








0 


Ct« 








Q C5 








w 








W 


Q 03 








pq 


0 


03 






pq 


S » 








CO 








CO 


E pq 








CO 


kO ro W 






CO 


H 








a: 


X 






oi 


H 


hd 






03 




M 






03 


H-J 


% 






w 


% 






w 


iJ ^ 


% 






pq 




V- ^ 






W 


h4 


rs 


0 


» 


> 


»-o 0 


W 


> 


hJ 


E 


0 


pq 


> 


hO 




0 


W 


> 


&-I 




r- 


3 


2S 




T— 


3 


25 


<! (Xi 




T— 


E 


z 


-caj 


Pm 




E 


z 


CJ M 


X 


T— 


2 : 


M 


X 


T— 


25 


H 


U H 


X 


T— 


25 


M 


CJ 


M 


T— 


Z 


M 








M 








H 










M 










H 








0 


Eh 






0 


Eh 








0 


EH 








0 


£H 








EH 


2 : 






EH 










EH 


z 








EH 


z 








0 


0 






0 


0 








0 


0 








0 


0 








0 


U 






0 


U 








0 


U 








0 


U 





M 

O EH 

EH 2S 

O O 
O U 



U 



O O 

^UCJCJ LOUUU 






o 

VO u u u 



o 

r- u u u 



o 

CO 



104 



CALL LAPLCE(X ( J, K) ,MU(MTYPE(J) ) ,STDV (HTYPE ( J) ) 
3,FX (K> ,IER) 



o 






























Eh 






















tsi 








CJ 


Csl 




























< 






2 
























Cm \ 




O 






























M 


























04 




H 






























CJ 
















n 








O 




2 
























t 






2 
















CO 








M CN 


O 




















04 








CU 


» 












2 








X 










CN 












O 








EH 






2 








C 


X 






H 








2 






O 




CD 




M 


2 






Eh 














M 




O 




Eh 








U 








> 






Eh 




4 




2 


\ 






2 








a 






U 


r5 n 






W 








ZD 








Eh 






2 








2 


04 






Cm 








CO 






2 


CO CO 






o 






















(m 


04 0L4 


Eh 




04 


o 














4* 








X X 


O 




X 


o 




















U 


Eh Eh 


<s 




w 


►4 






2 














H 


2 3C 


Cm 




1 








2 






o 








EH 








M 








O 








rn 






CO 


> ZD 


4* 










2 














H 


Q 2 














CD 








CO 






CD 


£H ^ 






ZD 








O 




2 


o 


04 






O 


CO 2 






o 












CO 


Eh 


X 






4 


2 






2 








1 




H 


O 


EH 






1 


^I. •< 






2 


w 






M 


04 


% 


CD 


2 






H 


Eh 


W 




W 


04 






.4 




NI 




•W' 






4 




04 




CQ 


X 






►4 


X 






ZD 






4 


O -V 


X 






H 






O 




2 


O 


2 






2 


t 


PH 




Q 


SD 






o 




04 










O 


o 


2 




W 


w’ 






2 


X 




t 


04 






2 


t 






X 








03 


2 


CO 


CO 


X 






2 


Eh f- 


2 




M 


SD 






W 




H 


2 


CO 






W 


« ^ 


2 




s: 








CO 


1 


2 


• 








CO 


O 








II 










2 




It 








CO « 


II 




w 






CO 


04 


Q 


2 








CO 








CO 








CO 




2 


W 








CO 


II X 






p:: 


bc3 






03 






M 








03 


2 






W 








w 


II 


.4 










CO 


Eh 


% 




> 


r5 o 


w 


> 




4 




n 


o 


CO 


> 


U II 


n 


CO 


2 






ZD 


2 


2 


*c 


Cm 




»— 


ZD 


2 


*c 




2 


M 


X 


T— 


2 


M 


04 


U 


M 


X 




2 


H 


Cm C4 


X 


2 








M 














M 








M 






o 


EH 












o 


Eh 








Eh 






EH 


2 












Eh 


2 








2 






O 


O 












O 


o 








o 




' 


o 


o 












CD 


u 








u 



w 



H W W 
EH O 
25 2 2 S 

o M H 03 

U EH Eh ZD 

2 2S Eh Q 

O O W 2 : 

O (J » W 



O 

U U U <7N U U U 



O 

o 

f- u u u 



o o o 

«- CN m 



o 



105 





4 ^ 


<U- 


4 !j> 






43 - 


'ill* 


4 -’,> 


43. 


4 i*-> 




a* 


4 ?^ 


a 


a 


a 


a 


a 


4 r;> 


a 


a 


a 


a 


a 


a 


a 


a 














a> 


■& 


< 2 * 


^ 5 * 


•i> 






a> 




a* 


3 * 


«£> 




a 


a 


a 




a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 




































































a 




































































a 














C’ 






































% 












0 




a 














1? 












X 


























cn 












in 




a 






















W 




0 






















X 
















z 




a 














•S' 








X 












Eh 














a 




II 












0 




a 




















Z 


Eh 




w 








X 














u 
















M 




a 














^1* 




2S 


0 






X 






X 


X 


















X 












EH 




a 


















0 


M 


in 




X 






z 


eh 














X 




X 
















a 


















H 


£H 


z 










EH 


X 














< 




X 












X 




a 




















Z 


X 




>* 








0 


















EH 












X 




ii- 














C' 




CJ 


CQ 


X 




X 






X 
















X 




H 












bH 




? 


















Z 


M 


EH 










0 


z 














X 
















M 






















0 


CC 


W 




X 








< 


















X 
















a 




















gH 


X 




z 






X 
















z 




0 












0 




<;• 




















C/^ 






0 






X 


X 






0 








X 




X 












0 




a 






z 








iJ 




Z 


M 


% 




M 






X 


M 






z 








> 
















0 




a 


















0 


Q 


X 




Eh 






X 








0 








M 




















■Vi 


















M 








U 






-a 


X 






M 








CJ 




0 






1 0 




























EH 


W 


w 




z 






> 








Eh 


















0 






Z 




a 


















0 


> 


> 




X 








«k 






X 








X 




X 






EH 






M 




i 








\ 






«0 




PC) 


M 


M 




X 






«c*; 


rn 






X 












11 












Z 




a 








0 






u 




M 


EH 


Eh 










X 








M 








z 




V 






X 






EH 




a 








CN 






iJ 




ca 


-SJ 


< 




z 






X 


II 






X 








M 




X 






-a: 






M 




a 








1 






ih 




6H 


h 4 


Z 




0 






M 








EH 


















X 






3 




a 








X 






43 . 




in 


3 


X 




M 






CJ 


X 






X 








X 










0 










a 








0 






4 SV 




M 


Z 


w 




Eh 






X 


X 






H 




X 


X 


X 




X 






X 






X 












0 






4 & 




Q 


s 


Eh 


t 


X 








X 






Q 








X 




Eh 










X 


CJ 




a 




X 
















CJ 


»-2 


5 m 


X 








EH 










z 


z 


X 




M 






X 




X 


X 




a 




X 




V 






05 «?• 




Z 




< 


EH 


H 






CJ 


M 






Z 




X 


X 


< 




X 






Q 




5 H 


X 




a 




5 M 


X 


X 






W ^ 




0 


ca 




H 


X 






z 








0 




> 


> 


> 




M 






0 




Eh 


> 




a 




Eh 


X 


bQ 






M <J 




m 


p 


X 




EH 






H 


X 






X 




M 


HI 






X 






U 






z 




a 














^•£1 




lA 




0 


H 


X 






Eh 


0 






X 




CJ 


CJ 


X 




< 


0 








X 


0 




a 














X 




H 


>H 




X M 






Z 


X 






M 








Q 




X 


X 




z 




0 


CJ 




a 












z 


^ i 3 * 




0 


EH 


w 


<J 


Q 






X 








0 




X 


X 


u 




0 


X 




X 










a 












1 


% 42 > 




cu 


M 


0 


X 








X 








X 




X 


X 






X 






X 




z 


Eh 




a 














W * 






in 


X 


0 


Q 






X 


t 














X 




X 


X 




EH 




0 


0 




a 












X 


Qj .* 4 . 




w 


z 


«: 


X 


Z 






X 


z 






Q 




z 


z 


X 






z 




X 




M 


z 




a 












X 


>H 4 , 1 ^ 




cr> 


w 


> 


X 


<x: 






X 


0 






X 




M 


M 


X 




X 


<j: 




X 




EH 






a 












X 






0 


Q 












X 


M 






X 








X 




0 


X 






Z 


a 


Q 




a 


















CXi 




<c 


X 


>H 






X 


EH 






H 




X 


X 


> 






-a: 




X 


X 


X 


HI 




a 












tl 






ca 


W 


CO 


X 6h 








X 






X 




X 


Q 


z 




X 


X 




X 


X 


X 


a 




a 














r: 1^- 




0 


X 


cn 


EH H 






X 


X 






X 




X 


CJ 


M 




X 








EH 


X 






a 






% 






X 


^«* 




04 




M 




X 






X 


H 






a 












Eh 


EH 




X 


X 


in 


z 




a 




X 


X 






Q 


X ^ 








U 


Z 


z 


t 




Eh 


X 










X 


X 


X 




X 


X 




z 


X 




X 




a 




X 


Q 






X 








a 


in 


w 


w 


X 




X 


Eh 






X 




z 


z 


z 




Z 


X 




-a: 




Q 


Eh 




a 




M 


X 


X 










<! 


M 


X 


> 


X 


z 




z 


X 






X 




X 


X 


X 






z 




X 


X 


M 


HI 




a 














z « 




ca 


X 


-a 


M 




0 




c 


M 


0 




Eh 




X 


X 


X 




X 


M 




-a: 


-a: 


X 


X 




a 












0 


Cn 43 . 




w 


EH 




CJ 


z 


M 




X 


Q 


X 








EH 


Eh 


eh 




< 






X 


z 


< 


0 




a 












0 


c/^ <: 




z 


H 


w 




0 


eh 




«: 




X 




X 




X 


X 


X 




X 


z 






X 


> 


CJ 




a 












0 


H «. 




w 


W 


X 


X 


X 


< 




X 


Z 


H 




0 




X 


X 


X 






-a: 




X 


0 


z 


X 




a 












H- 


O ^1* 




CJ 




Eh 


X 


X 


X 






0 


X 
















EH 






0 


z 


H 


< 




a 














0^ «3* 






cn 




X 


H 


0 




EH 


X 


Z 




Z 












X 


X 




X 










a 














< 3 * 






z 


Z 


-a; 0 


X 




X 


X 


< 




< 




f— 


CN 


m 




X 


H 




X 


1 


1 


1 




a 












0 


CO * 






X 


w 


> 


X 






X 


M 


X 




X 












Eh 






X 










a 












Eh 


Z i2* 




in 


0 


> 






z 




z 


0 






z 




M 








X 


X 






0 




m 




a 


cr 








0 


X 


M ^ 




M 


EH 


M 


-a 


X 


0 




M 


X 


X 
















0 


X 




1 




1 


1 




a 


a 








f 


0 








W 


CJ 


X 


w 


M 












I 




X 
























a 


X 




% 




0 




3 




Z 


Z 




X 


EH 


X 




f 












X 








1 






X 










a 


X 




X 






ZD 


O iS* 




cn 




W 


H 


X 


X 












X 




X 








X 






X 










a 


CJ 




Q 


% 


II 


Z 


« <-3* 




cn 


Eh 


X 


u 


X 


X 




X 








z 




Eh 








X 






M 










a 


X 


M X 


CJ 


< 






CQ ^ 




w 


-a 


X 


X 


z 


u 








































a 


EH 


«5 






X 




o 3 




0 


X 


< 


X 0 


X 








































-■'% 


z 


X 






Q 


X 


cn 3. 




04 


EH 


> 


f< 


CJ 


X 








































a 


M 


X 






X 


M 


-t? 






















































a 




a 


a 


a 






14 






















































w 














■"S' 






















































a 
















43* 


*& 








4:> 


a* 




‘3' 


e* 




a* 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


a 


c 


<5 


a 


























•a* 


a* 


w 


a* 


a- 


-Pf 


a 


a 


a 


a 


a 


a 




a 


a 


a 


a 


a" 


*i5 


a 


a 


a 














CJ 


CJ 


u 


CJ 


u 


u 


CJ 


CJ 


CJ 


CJ 


u 


CJ 


CJ 


u 


u 


u 


u 


CJ 


CJ 


CJ 


u 


CJ 


CJ 


CJ 


CJ 


u 


CJ 


CJ 















106 



c 



a 



o 

o 

o 

















O 






04 






. 












cr 






























43* 








O 






O 




O 








O 






=3- 






cr 










n 


t 


o 


O 












o 






y? 


o 


vD 


VO 


o 






O 




CD 




















Eh 










!l 


II 


O 


o 


o 






O 














H 


£-i CO 








CD 




X 






(t, 




O 


o 04 












ii4 






Q 


a 


O 


C) X 




EL4 






Cu 








0( 


04 






X 


Q 


<T> 


X 


Q 0> 


t 
















04 


o> 




u o> 


Eh 












m ^ 


« 






f 




hO 


c 










o 




11 


O 


Eh 


M O 


t 




o 


5^ 


CO 


t 


t m 


OJ 




Eh 




Eh 




ii 


EH 






H 


Eh % 


« 


X 


O 


t 


X O 


Cu 




O 


• 


• 


OJ 


CD O 




Cui 


CD 




Cu CD 


Q 


X 


CD 


O E-^ 


W 


• 


# <N 


c 










a 






t CD 


hO 






















9 


t 



Cli 



04 







r* 


w 


W 


O 


W 




PO 




W 


w 








Q 






















04 


04 


T— 


HD 


Ct4 


o 




3 0^ 


o 




Cui 


04 


U 
















CX4 




T— 


X 


X 






M 


25 


H 


2 M 


2 




a 


O 


















Q 




II 


H 


Eh 




M 




W 




M 


M 


C 


04 


04 


II 


w 


rn 




W 


f— 


CO 




04 


o 


M 




N«.^ 


o 


£h 




Eh 






Eh 




N»l» 


'W 




=D 


1 




2 


1 


o> 






t 








Eh 


25 




CS 




2 


2 


II 






Ct4 


2 






2 




2 


2 


11 


o 


O 


Elh 


Ct4 


O 


o 




O 




O 


o 




CX4 


Cl. 


a 


M 


It 




M 


II 


H 


03 






in 


H 


H 


CD 


o 




O 




u 


u 


< 


H 


H 


CJ 






O 








3 


Pu 


II 






























2 


cr; 


Eh 


2 


C35 


2 




a 




o 




























O 


W 


O 


o 


W 


O 


CO 


CJ 


< 


o 




























CJ 


H 


CD 


CJ 


H CJ 


03 














o 




O 




o 


o 










o 






o 




CT> 
















r-» 




(N 




rn 












m 






VO 









107 



END 



SUBROUTINE NEGBIN ( X, M U ,S TD V , T YPE ,FX , lER) 



& 


‘a' «ji 


Ji' 


Ci> 




c> 


•S' 


•S 


<C> i'f ', a Hr 


'S 




<s 




^ ‘Cl> 


•s 


<r? 


«s 


«(ii ‘Si» 




•S 


fl- 






J 






c> 


•S3* 


'S 


&•»»■£!■ 


-s 


•s 


^3* 




<» o 


•» 






<3> w 


O -a 




'S 







>i 




z 






















% 








0 

cn 






1-1 




z 




















Z 


m 








z 






W 




o 




















Q 










O 


ss 




> 




H 




z 


cn 














U 


II 








H 


O 




M 




Eh 




X 


M 
























EH 


H 




H 




CJ 




Eh 
















Z 


z 








-*S 






<C 




z 






X 














-C 


z 








z 


O 


w 


Z 




z 




Z 


















X 








z 


CQ 


> 


a 




z 




o 


% 






z 








X 










EH 


M 


H 


u 










no 






-a; 








z 










M 


« 


H 






z 




Z 






0 


M 










z 










H 


•eC 






o 




z 


II 




z 


X 








z 


O 








Z 


cn 




<r: 


t 


H 




z 






o 


o 








z 


z 








o 


H 


2 ) 




>* 


EH 




z 


Z 




H 


z 








> 










o 


a 


z: 


cn 


EH 


Z 




rtj 


z 




Eh 


M 








M 










t* 




3 


O 


H 


Z 




> 


X 




Z 


z 








iD 


0 


0 # 










u 




Z 


H 






Eh 




z 














o 






z 


c 




% 


M 


Z 




<c 






H 


z 








X 


X 


Eh 






H 


M 


PC 3 


u 


Z 


Eh 




in 


Z 




Z 


> 










II 








X 


z: 


O 


z 


-c 


cn 




in 


O 




EH 


H 








Z 


V 


Z 






Eh 


o 




h 4 


z 


M 




H 


z 




Z 


Eh 








H 


X 


-c 






M 


:z 


>* 


<C 


o 


Q 




u 






M 


< 












z 






Z 


M 


e-i 


> 


z 






in 






Q 


O 








Z 




o 








CQ 


M 




z 


Q 




z 


• 






Z 




X 


X 


Z 


X 


z 






z 




CT} 


<< 




z 




<i: 


z 




Z 


Z 








z 


Eh 






z 


iD 


u 


Z 


cn 


w 


<a: 






o 




< 






z 


Z 


z 


H 


z 




z 


Z 


> 


W 


cn 


z 






o 


H 




M 


z 




z 


Z 


-c 


Z 


G 




X 


z 


H 


Q 


H 


Eh 


>H 




z 


H 




C 


X 




> 


> 


> 


M 


o 




Eh 


> 






U 




£H 




H 


z 




o 


Eh 




M 


H 




Z 


u 






z 


-S 3 


W 


cn 


z 


M 




Eh 


z 




Z 






O 


O 


z 


C 0 






Z 


o 


O 


Z 


z 


z 


in 




Z 


M 




H 


z 








Q 


X z 


z 




O 


u 


W 


E-* 


<i: 


> 


z 


t 


Z 


Z 




Z 


o 




X 


X 


u 


o w 


z 








23 






M 


w 


cn 


cn 


Eh 










z 


z 




Z Eh 


z 




z 


ch 




cn 


w 


O 


z 


z 


z 


Z 




z 


z 








z 


Z Z 


EH 




o 


o 


U 


w 


z 






o 


z 


M 


• 


> 


o 




z 


z 


z 


X 


z 




M 


z 


CO 


s 


Eh 


W 


h 4 


H 


Oi 


Q 


Z 


H 


M 




M 


M 


z 


z -n 


z 




Eh 




O 






z 


< 


H 


z 




z 


Eh 


Eh 








z 


O z 




z 


U 


G 


Oi 


M 


Z 




H 


<t: 


z 


Z 


Eh 


-x: 


-s: 




z 


Z 


> 


c 


« 


z 


Z 


M 


cn 


w 


u 


-c 


Z 


z 




-C 


Z 


z 


H 


t 


a 


Q 


z 


z z 


z 


z 


Z 


Q 


Z) 




> 


> 


o 


o 


z 


M 


z 


z 


> 


z 


z 


U 


H 


z 


EH 


£h 


z 




Oj 


cn 


M 




z 


w 


z 


S 3 




z 


Z 


o 








Eh Eh 


Z 


Z 


in 


X 




z 


CJ) 


<c 


H 




Eh 


O 


z 




Q 


H 


z 


Z 


Z 


z z 


X 


Z 




X 


1 ^ 


cn 




cn 


Z 


z 


Z 


z 


c 


z 




EH 


z 


Z 


z 


X z 


-c 




Q 


Eh 


<C 


z 


W 


z 




o 


S 


M 


z 


X 


Q 


Z 


z 


z 


z 


<C z 


z 


Z 


M 


H 


C £3 


c-i 


z 


M 


w 


M 


<C 


Z 




Eh 


Z 


z 


z 


z 


z 


Z M 




-C 


z 


Z 


W 


w 




u 


> 


cn 


Z 




H 




<d 


H 


E-- 


EH 


E-- 


c 


z 


X 




O 


z 


cn 


< 


in 


M 


z 


-C 


Z 


z 


z 


Q 


Z 


Z 


Z 


Z 


z z 




z 


> 


o 


w 




> 


z 


Eh 


z 


Z 


> 


z 


o 


Z 


EH 


z 


Z 


z 


-s: 


z 


o 


z 


z 


o 






•ci: 


< 


CJ 




M 


Eh 




< 


Z 








Eh 


o 


z 


M 


c 




-s; 


z 




Z 


z 


H 


Eh 


Z 


z 


EH 


H 








z cn 


z 








<r: 


z 


o 


w 


w 


z 


Z 


Z 


o 


-c 


Z 


Z 




CN 


m 


Z M 


z 


1 


1 


1 




E- 


H 


z 


z 




Z 


iD 




z 












Eh 


z 








U1 






EH 




z 


z 


Z 


z 


X 


1 




11 






Z X 




o 


t* 


ro 


M 


Z 


z 




cn 


o 


M 


Z 


<c 














o z 


1 




1 


1 




o 


z 


cn 


w 










1 


> 




z 
















z 


M 


M 


z 


£H 


w 


1 








Q 




z 






1 


z 








M 




Z 


z 


z 


cn 








z 


H 




X 






X 


z 








CQ 


u 


6 -« 


z 


z 


z 


X 






X 


Z 




Eh 






z 


w 








a 


Z 


cn 


E-i 


z 
































w 


z 


M 


W 


o 


>H 






























z 


Cm 


Q 


z 


u 


z 





































•s 




•el> 


-s 


•& 


<s 


•s 




C 




-s 


•s 


•s 




e 


•S 


•& 


•S 


4-> 




-s 


e 




•s 




« 












15^ 






•£3* 


<3* 




VJ’ 




•£;* 




<3* 






w 


<*• 




<3 




i> 








U 


u 


u 


u 


CJ 


u 


u 


CJ 


CJ 


CJ 


CJ 


u 


U 


u 


u 


o 


U 


U 


u 


CJ 


u 


U 


U 


u 


u 


U 


u 


u 



103 



2?^ l^. * S${ 2^4 I'y SgK Ijl 2s JCi 5^ 5^ C - ^ C ■ C* ^ JSJ ■? •> «0^ 

INTEGER^S^U 

I, lER, TYPE 

REAL‘S 4 

$ CDF, PDF, FX, HU, 






0 

CN 

1 

w 

o 



V 

% 

X (SI 



Pm 

0 

01 



0 

£ 

1 

w 

S 

o 

<2* 



cc 

< 

05 > 

(N 

i- > t *5 

> N 



CN O Cm 



o 

o 

o 



f f ro W 
Eh Eh % C3 Pm 
O O Z M 
• • CS M 

^ &H 

- w pLl o s 
H O 4 p4 r- O 

- X X CJ 



w 

25 

M 

EH 

2: 

O 

U 



Pm 

H 



W 

2: 

H 

£H 

2; 

O 

U 



M 

25 





























0 


% ^ 




























«< U 5 








0 






0 






0 






1 










cr 












Eh 




























0 
















0 






0 






0 
















Eh 






EH 






















0 






0 












0 










0 






0 






X 






< • 




0 


0 
















Pm 






^0 


\ 


VO 


vO 






Pm 






O4 












0 








X 


Q 


a\ 


X 


Q 


C 7 > 


f 




o> 


II 11 


f 


0 


0 






04 


cr> 




U 


(T> 


Eh 




(T* 




r— 


Eh 


Eh 


w 


• 






» 






»-3 


<3 




Pm Cm 


X ^ 


0 


0 


04 


Eh 


)l 


0 


EH 


)l 


0 


f 




0 


Q Q 


05 &3 3 


0 


0 


X 






Eh 


1-3 




EH 




II 


Eh 


04 04 


«: 2; S 






Eh 


• 


X 


0 


f 


X 


0 






0 




> 0 






% 




Cm 


e) 




Cm 


0 


0 


X 


0 




1 


r- 


cn 




< 












u 






h 4 








0 




















S 3 (S 3 



Cm 

Q 

CI4 



W 

2: 

O 



w 

1-3 

f 



pq 

=) 

25 



Eh 
O 

Cm 
Q 

Cm CJ 

a a 



% 


Eh 


0 C 


2: 


0 




11 


EH Eh 




H 




04 


Oi 


11 


cy 


rn 


> 


% cn 


S 3 ^ 


03 


O4 


0 


M 




0 


Eh 










0 


1 


Q 


0 


0 






» 






6 -* 


2: 


II 






Cm 


2: 




EH 


az II 


II E 


II 


II 


0 


0 


Pm Cm 


0 


0 




Cm 


O4 


Q 


M 


)l 


cn 


03 










m 


H H 


CD 


U 




M 


H 


U 


EH 






03 


0 I) 


Pm 


O4 


II 


















2; 


03 




< 


K 


0 


Q 




0 
















0 


w 




> 


03 U 3 


O4 


U 


< 


Q 
















a 


M 






0 


0 


0 


0 


0 




CN 


rn 


cr 


uo 



109 



GOTO 99 



CONTINUE 
lER = -1 
CONTINUE 
RETURN 
END 



O 

vX> 



110 




C0H41NaB 



SUBROUTINE LOGIST ( X, MU ,S TD V ,T YPE ,FX , I ER) 



■& a- t. 




«■ 




^r> 




i? 


•S' 


‘J3- 


'S’ 


i'> 


«C2- 


•{S' <I> 


O il ■£> 


«• 


<3» «& 










•s* 




{J 


{} 


•d* 


^S* o 

















X 












W 










EH 












X 










Z 


X 








z 


£H 










X 


o 






o 


o 






X 






z 








H 


M 


CO 




o 


z 




o 


z 






6-* 


Eh 


z 










z 


o 






U 


ED 


X 




X 








M 






z 


PQ 


X 




X 


X 




X 


£H 






3 


M 


EH 




X 


H 




a 


X 






Cu 


03 


W 












X 








H 


X 




> 


X 




X 


M 






z 


V) 












o 


X 






o 


M 






*< 


% 




X 


EH 






H 


Q 


>H 




X 


rn 




£H 


X 






H 








X 






X 


H 






3 


W 


M 




H 


It 




z 


X 






PQ 


> 


> 




U 














M 


H 


M 




X 


X 




X 


u 






ce: 


Eh 


EH 




X 


X 




o 


H 






H 


•< 


«< 




*< 


X 






EH 


t 




cn 


hJ 


z 






Eh 




z 


X 


X 




H 


3 


X 




X 






o 


M 


eh 




O 


ac 


w 




X 


X 




M 


CD 


<c 


X 




3 


Eh 


« 


Eh 


O 




Eh 


O 


z 




O 


U 




>» 




X 




X 


X 


M 


z 


H 




-a: 


Eh 


CD 






X 




EH 


X 


Bri 


0:5 




M 


Z 






H 


X 


X 


> 


in 


O 


X 




H 


t 




X 


X 


X 


M 


M 




O 


M 




z 




EH 


EH 




CD 


O 


>H 




X 


Z 


O 




X 




X 




o 


Eh 


W 


< 


X 


H 




H 


X 


X 


X 




H 


X 


X 


X 


Eh 




Q 


o 


z 


X 




CO 


hJ 


O 


X 


X 












Ce3 


Z 


-a: 


X 


X 


X 




CJ 




X 


z 


tn 


w 


> 


X 


X 


M 




M 


o 


o 


H 


o 


Q 






X 


X 




EH 


M 






Oj 






X 


X 


EH 


• 


X 


eh 


X 


X 


05 


W 


in 


X 




X 


X 


H f 


*< 


Eh 


X 


z> 


X 


CO 


eh 


X 


M 


X 


CD X 


M 


z 


X 


04 


Eh 


M 




X 


a 


s 


O Eh 


> 


X 








U 


z 


eh 




X ^ 


X 


z 


X 


hJ 


X 


CO 


X 


X 


u 


z 


Z 


Q 


o 


z 






X 


> 


X 


H 


<c 


X H 




z 


X 


CC5 


X 


•< 


M 


-a: 


Eh 


X 


X EH 


X 




X 


W 


H 




CD 


X 


X 


rf: 


Eh X 


X 


X 


EH 


Z 


M 


w 




-a; 


H 


X 


X 


rf: 


o 


X 


w 


W 


X 


X 


X 


CD 




X 


X 




X 


a 




Eh 


X 




O 


EH 


O X 


z 


X 






CO 




X 


Eh 


X 


X 


X 


rt: 


c 




c 


z 


Z 


<; 


X 




X 


z z 


Eh 


X 


r-* 




X 


w 


> 


X 


X 


Eh 


c 


X 


Eh 




tn 


X 


> 




z 


X 


o 


X X 




X 


it 


H 


£H 


H 


< 


H 


EH 


o 


z o 


I 


Z 






W 


CD 


in 










> 




X 


EH 


X 




in 


1 






1 


X 




X 


CO 




W 


H 








X 


Eh 




X 


M 


EH 


X 


U 


X 






z 


X 






O 


C 


1-3 


CO 
















O 


X 


-a: 


X 
















►-? 


EH 


> 


«*: 



















•S' 


•d* 










•o 




•d' 


•G 






lii’ 










‘d' 














•2^ 


'd- 






•S' 




•d' 


u u 


u 


U 


u 


u 


u 


u 


CJ 


CJ 


u 


u 


u 


u 


u 


u 


U 



<SV 




4Tv 


i? 






•ili» 


















43> 


^ i2> 


































•?*. 


























«> 










































% 
























X 


m 










•O 














X 












‘d’ 














u 


tt 










<3* 


































% 






X 


X 
























< 


X 


















o 








X 


















• 






X 


H 










<3- 








t— 






X 












■0 








\ 








X 










•(S' 








X 






z 


o 










4S- 








z 






X 


X 










<!’ 








c 


X 




> 












e 














M 


























CD 


t 




f « 






•vi 




















o 




















X 


X 




EH 






















It 










^3* 








% 






z 


V 




X 














\ 






M 


X 




*< 






<£> 








m 












X 






•S' 








t 






X 






o 














o 


% 


X 


X 


X 




X 














\ 


M 




X 


EH 








X 










X 


CU 


z 


X 


H 




X 




X 










X 


X 


X 


c 


X 




Q 




X 










c 


o 


> 


> 


M 




O 




EH 










X 


CO 


M 




X 




CJ 






i. 












CD 


X 


c 


* 






X 
















X 


X 


X 


z 




o 














X 


CJ 


o 


X 


X 






i. 












X 




X 


EH 


X 




z 


•C^ 














X 


X 


X 


EH 




o 














z 


X 




z 


X 




M 














H 


X 


X 


c 


X 




EH 
















X 


o 


X 




z 


CJ 














X 


> 




c 


X 


X 


X 


'w' 












X 


z 


X 


X 


X 


X 


X 












X 


CJ 


M 


X 




£H 


Eh 


X 


i. 




X 






z > 






X 


EH 


X 


X 


CO 


♦S. 




CU 




% 


CJ Q 


X 


X 


X 


X 


z 


X 








X 




X 


X EH 


z 


z 


z 


X 


c 




a 


C 




Eh 




X 


CO CO 


X 


X 


c 


z 


X 


X 


H 














X 


X 


X 


w 




c 


X 














EH 


Eh 


c 




X 


z 


-a: 












V 


X 


X 


X 


z 




X 


> 












\D 


X 


X 




-a: 


X 


o 


z 












CN 






EH 




o 


z 


H 












cn 






X 


CO 


X 
















in 


CM 


m 


X 


M 


X 


I 


I 


















H 




X 






-d 
















X 


X 




o 




& 










r— 






o 


X 


t 




t 


? 










» 
















S’ 


X 


% 






rn 






t 




X 






tJ* 


X 


X 




CD 


\ 






X 




X 






* 


X 


X 


S’ 


X 


w o 






X 




H 






d* 


X 


H 


X 


-a; 


X z 
















-S’ 


Eh 




•< 




















S 


Z 




X 




















S- 


M 




X 







^ ^ ^ 



«• {J ^ i> il- <•> '5!' {3- 

^ lii* w <5« i3* 

uuuuuuuuuuo 



111 



INITIALIZE CONSTANTS 



s: 

Q 3 

O 

H 

in cu 

z 

< X 

(X 

en ^ 

























3 


X 






























< 


3 


















> 












3 




















Q 












CD 


1 


















H 












W 




















CO 












Eh 


w 






























z 


z 








A 






















M 


o 








> 
































a 










m 












X 










H 










Eh 












EH 


X 








CQ 










03 












H 




















O 












h4 


X 


















CO 












M 


3 






























3 










m 


















CD 




<C 




























3 




3 


o 








(X 










o 








»c2: 




O 


h4 








o 








CX 


t 






fX 






3 










cn 








Q 


— cr 






Q 


3 




CU 


i2> 
















CX 


O 






u 


Z 






















cc 








< 




3 


m 








X 








Z 


X 






Z 


Eh 




O 


Eh 
















o 








O 






3 


3 








A 








M 


3 <N 






H 


4* 




Pu 


a 








=5 








PH 


CO i% 






EH 








CO 








s 


o 


O 




3 


O ^ 






3 


W 




CO 












cr 


cr 




tri 


O X 






3 


z 




w 










1 








M 


3 






M 


o 




3 












o 


o w 




cc 


X U 






3 






1-4 


t> 








X 




EH CU 




Eh 


w 






Eh 






< 


Q 










o 


O » 




cn 


W CO 






CO 






> 


Eh 










o 


O H 




M 


z 






M 






CO 
















Q 


o ^ 






Q 






X 










o 




% 














3 






H- 






r— 


# H 


r— 


on — . 






II M 


o 




3 


<< o 




W 




O 




1 


ro o< 




o 




Eh 


CU 


cn 




Eh 


3 to 




Eh 


3 


m 








• 


• m 




3 


X 






3 






3 


3 






II 


H 'iS- 




e ^ 




Ch 


3 II 


O 




CU 


II O 




3 




o 






X 


h-3 


o O 


W 


S 


O 


Eh 


» 


3 


EH 


w 


fi 


II 


£H 


W 


3 


O* {X 


f 


• CN 


3 


o 


W X 


O 


3 


O 


X o 


3 


O 




O 


3 


W 


cn 1-3 




% 


Z 


u 


CO 


C5 


Z 


U 


C*H CD 


Z 


U 


X 


CD 


Z 


X 


«d 


W 


w o 


M 








H 






M 








M 




li Z 


Oi 


CU t- 


Eh 








E-» 






£-1 








Eh 






X 


>» 


Z 








Z 






Z 








Z 




m II 


Eh 


Eh 


O 








O 






O 








O 










U 








o 






u 








CJ 




X O 




Eh 




























O « 


PX 


Cl, O 




























cn -a: 


M 


H O 





























w 



z 

o 

u 



u o u 



o 

f- o u u 



o o 

CNOUO OOUOU 



O O 
^ to 



112 



RETURN 

END 





4r? 






•I3* 






'C? 






e 


"tl* 




<£: 




‘Si’ 


42> 


<12* 






•& 


4ii> 




& 


«& 


<& 


‘O;* 


•» 












O 








fj* 






{!« 


<2- 










«2- 




‘Si* 






■& 


















‘S' 








































































































































































































c 










































% 






























W 














3 
















3 




m 














e 
















3 














Eh 
















O 
































2 
















3 






3 










a 




It 


























z 


O 








3 








3 






O 








































o 


H 


CO 






O 








S 
















3 




3 


























H 


H 


3 








3 






o 






z 










rt: 




3 




























3 


3 






W 


»cd 






3 






o 














X 














^ir 












U 


3 


3 






3 














H 










X 




EH 














it 












z 


M 


EH 








3 






3 






EH 










3 






























3 


3 


3 






rt: 


H 






O 






3 














3 


























b 




3 






> 














3 










z 




O 














‘frS 








<:• 






CO 










X 






Q 






H 










3 




3 






















c 




Z 


M 


% 






rt: 








O 






3 










> 


























vi 




o 


Q 


>4 






3 








3 






EH 










H 


















.fj 












H 




»4 






3 


m 






EH 






3 










O 




t 






« • 




















M 


W 


W 






M 








3 






H 




















o 




















3 


> 


> 






CJ 


It 






s 






Q 










X 




X 






EH 








‘i2’ 








t 




3 


H 


H 






3 




























II 






















'kt 




H 


H 


EH 






3 


3 






3 






3 










z 




V 






3 




















3 


< 








< 


3 






O 






U 


t 








M 




X 






< 
















03 






3 


3 








>4 












•a: 


3 


















Z 
















w .» 




CO 


3 


3 






W 


EH 






3 






3 


EH 








3 










Oi 
















H 




H 


S 


W 






3 








o 






3 


C 




X 


X 


3 




X 






3 








4n» 












Q 


3 


Eh 


• 




EH 


3 






H 






rt: 


C 








Z 




EH 










3 


















o 


•-0 


>4 






O 






Eh 






3 


H 




z 


z 


3 




M 






3 




3 












Pm -w- 




W 




lo: 


Eh 




O 


3 






3 








£h 




3 


3 


kC 




3 






a 




X 




•2* 












U 


03 




H 




3 








3 






3 


3 




> 


> 


> 




H 






o 




Eh 












W 




c 


O 


3 






M 








H 






3 


3 




H 


H 






3 






u 








V* 








CU -?} 








C 


H 






• 






3 






EH 






O 


O 


3 




< 


t 








3 












X -ij 




04 


>4 




3 




3 


3 






EH 








3 








Q 




3 


3 




z 




O 












H i:> 




rt: 




w 


«C 




3 


O 






3 






3 


3 




X 


X 


CJ 




O 


3 




3 








•r: 












3 


H 


3 


3 




3 


H 






H 






O 


Z 




3 


3 






3 


Eh 




z 




z 




e 








> 4it 






CO 




O 




3 


H 






O 
















3 




3 


3 




EH 




O 












Q 




w 


2: 


< 


3 




3 


3 












z 


3 




Z 


Z 


3 






Z 




3 




M 




4% 








H =» 






W 


> 


3 




04 


3 






3 






O 


O 




H 


W 


3 




3 


C 




3 




EH 












C/) ^ 




o 


Q 








W 


H 






a 






M 










3 




O 


3 






z 


U 
















04 




•< 


W 




3 


3 


• 










Eh 


3 




3 


3 


> 






C 




3 


3 


3 




<1» 








3 ^ 




3 


3 


CO 


3 






EH 


3 




3 


• 




< 


EH 




Q 


Q 


Z 




3 


3 




3 


Z 


3 




*£}> 








S 




3 


3 


CO 


EH 




3 


3 


3 




3 


3 




H 


Z 




3 


U 


M 




3 






Eh 


£h 


3 




<2> 




3 








04 




M 






3 


H 


Eh 




C 


Eh 




> 


3 












Eh 


H 




3 


3 


CO 








3 












CJ 3 




Eh 


Q 


3 




3 


*< 




3 


Z 




3 


3 


3 




3 


Z 




Z 


3 






'IJ 




X 








3 




CO 


W 




3 




E 






3 




Q 


o 




Z 


Z 


Z 




Z 


3 








Q 




<s 




Eh 




fj* 




»cC 


U 


3 


> 




3 


3 


<< 




3 


H 






z 




3 


3 


3 




< 


Z 




3 


3 


H 




O' 








» «> 




03 


3 


<U 


M 






CJ 


3 




3 


EH 




Q 






3 


Z 


3 




3 


M 




4«32 


-a: 


3 




r 








U 




w 


Sh 




CD 




3 


f< 


-i: 




EH 


3 




3 


3 




Eh 


EH 


EH 










3 


z 


<d 












ti^ 




21 


H 


W 






-C 


3 


3 






3 




-a: 


O 




3 


3 


3 




3 


Z 






3 


> 




42 








Oi <f!2. 




CjO 


W 


3 


W 




3 


3 






3 






Q 






3 


3 


3 






rt: 




3 


O 


z 




42* 








^ e 




o 




EH 


3 








EH 




O 


3 




Z 


Q 












Eh 






O 


z 


M 




43* 








>-i ii- 






CO 




hO 




Eh 


3 


3 






3 




< 


Q 












Z 


3 




3 








«K> 








<T> 




< 


3 


3 






3 




3 




25 


3 




EH 


3 






CN 


no 




3 


H 




3 


1 


1 




<.1* 








w 2^ 






3 


W 


> 




3 


3 


H 




rt: 






3 


H 












Eh 






3 








‘£r 








!25 -S- 




CO 


3 


> 






3 


3 


3 




3 


3 






3 




It 








Z 


X 






o 






4^ 








H ^ 




H 


EH 


H 


<< 




H 




O 




S 


O 




1 


Z 












O 


3 




1 




1 












H 






W 


c:) 


3 


















> 






3 






















4-7 


3 


% 




3 <:• 




w 


3 




3 




1 








1 






Q 






3 








1 






3 








45* 


3 


3 




O <2 




u 




w 


H 












3 






EH 






X 








X 






3 








•S' 


CD 


3 


<2 


03 ^ 




3 


EH 


3 


O 




X 








3 






3 






£H 








3 






3 








•iil* 


3 


H 


3 






04 


0SC 


h4 


3 
















































Eh 




-i: 


3 




<< 


3 




3 














































i2* 


Z 




3 


cn •& 






EH 


> 


•< 
















































H 




3 






























































42* 




























































'w 








i- 
























































12* 










43* 






e- 


'S* 


e 


<4 


«& 


«• 


e 




4^ 


it 


v> 
















<2* 


<* 


c 


‘S' 




42* 












iS- 


'ft' 




«* 


* 




i§* 


iJ* 


<2Sf 


Q 




<2* 




i> 














4^ 






ik 


C' 


it 


<3* 










U 


a 


CJ 


O 


U 


a 


a 


U 


a 


CJ 


U 


U 


U 


u 


u 


U 


U 


a 


a 


CJ 


a 


a 


a 


CJ 


u 


u 


u 


CJ 


CJ 









113 



ARG, FX, ONE/1.0/, SQR 2/1 . 4 1 4 2 1 35/ 

TWO/2. 0/, X, ZERO/0.0/, MU, 

STDV 



INITIALIZE CONSTANTS 



CN 
« 
Of 
CO 
CN \ 



CJ u u 























z 


05 ^ 






















z 


Oi X 






















o 


CO Du 






















Du 
























CO 


\ ^ 






















z 
























< 
























05 


X 2B 






















EH 


tX4 EH 






















hO 


1 






















-a: 
























03 


O O 






















CJ 


Z Z 






















DO 


EH EH 






















EH 
























Z 


CJ CJ 






















M 


O O 
























hO hO 






















X 
























Eh 


^ w 






















H 




















X 




hO 


> > 










> 








04 




M 


Q Q 










Q 












CQ 


EH EH 










Eh 








1 




C 


in in 










in 












CQ 




















w 




O 


-f 1 








Cl< 






Du 




z 




Z 




> 






a 






Q 




o 




04 


ZD ZD 


a 






04 


CN 




U 










z z 


Eh 








05 








II 




z 




CO 






2 


O 




Z 


O 






o 


tl II 








o 


in 




o 


CS 


X 




« 




V 






M 






M 


EH 


Cu 




04 


X X 








Ph 






EH 














O 


















CO 




o cr 






00 






CQ 




o 




w 


in in 


s: 






M 






H 




03 




Z5 


t t 


O 


O 


W 


05 


CJ 




05 


CJ 


W 




hO 


o o 


1 Eh 


H 


04 


EH 


cn 




EH 


04 


IS3 




<s3 




O 


O 


X 


CO 


< 




CO 








> 


t t 


X O 


CJ 


H 


M 


1 




H 


1 


• 






EH PO 








Q 


w 




Q 




w 




X 


iJ CJ 




— 






04 






04 


CJ 






t t 


CQ r» 


no 




W 


X o 




W 


X 


• o 




DO 




-aj 




O 


EH 


CO m 




EH 


w 


in 




EH 


X X 


• 


t 


ro 


TD 






Z 




X 




ZD 


pL, 04 


^ EH 


Eh 




04 


H O 




04 


II 






04 


W 


H) 


CJ 


o w 


s: 


eh 


w 


z 




Eh 


W 


z 




CN t 


« 


nj r> 


o 


X o 


ZD 


o 


X 


Pt4 o 


Z) 


o 


04 04 


05 




% :zs 


u 


Du CJ 


2: 


u 


04 


Hi CJ 


z 


CJ 


H H 


o w 


W 


O M 






H 








M 






in 04 


04 


T- EH 






Eh 








Eh 






X 


>4 


Z 






Z 








Z 








EH 


O 






O 








O 










o u 






U 








U 






o 




Eh 




















05 Pu 


04 


O 




















H 


H 


CJ 
























O 






O 








o 












c 

c 


U 


CN 


c 

c 

c 






cn 


c 

c 


U 



o 

in 

O 

O 



I 
)l 

W 05 
W 
:z 
Hi 
Eh 
25 

O 
U 



o 

; 3 - 



W 

n 

:z z 
H C3 
EH S 
» Eh 
O W 
CJ 05 



O 

in 



114 



END 






i 







«d* 


<rj. 


G> 


«G» 


<£ 


42 


42 


4-i. 


42 


4;? 


42 


4-,*’ 


42 




42 


42 


e 


4Ti> 


42 


45, 


42 


4!i> 


42 




42 


42 


42 




42 




4:e 


42 


42 


42 




<2* 

•fiV 




<3* 






O' 


<2 


'2* 


4> 




42 


i2 


42 


42 


<2 


<2 


42 


42 


12 


42 


42 


42 


42 

p n 


42 


42 


42 


42 


42 


42 


42 


© 


42 


42 


42 


12 




ii' 

Ci- 












































Z 

H 
































X 




f 








a 




























Z 
















z 




Q 






cn 


•a* 




hJ 




CO 








z 












Z 
















Z 
















< 




Z 






s: 






a 




Eh 








Eh 












z 


Z 








Z 




f 


X 




Z 












z 




Eh 






% 






CO 


lO 










H 










C 


z 


Z 




O 


Z 


Z 






< 




X 












z 




U 






Eh 


•?> 




Eh 


o 


Z 


a 




O 


Q 










z 


c 


a 




z 


z 


Z 




Eh 


z 




Eh 












z 




Z 






C 






0^ 




CO 


w 




Eh 


Z 






a 




z 


H 


z 


Q 


z 


z 


z 




X 


z 
















z 




z 






3= 


45* 




CO 


CO 


w 


z 












z 




H 


z 


<< 


Z 


o 


z 




Z 


CJ 


z 




% 












£ 




X 


X 




Ch 






04 


z 


z 


H 




Z 


z 




Z 


Eh 


Z 


U 


c 


z 


< 


z 


z 


z 


Z 


H 






Z 


Z 














z 


z 




% 


'Si- 




CO 


o 




«c 




z 


z 




CJ 


ZD 


Z 


Z 


> 






s 


CJ 


z 


Eh 


Z 


z 


Z 


Z 


X 










z 






X 




X 


*» 




CC3 


M 


z 


Eh 


X 


M 


Eh 




< 


Z 


z 


Z 




o 


z 


o 


CJ 


z 






X 


X 


M 


EH 














H 


z 




pL4 


•S2* 






Eh 


z 


Z 


Z 


Eh 






z 


z 


z 


<< 


z 


EH 


o 


u 


ZD 


z 


Z 


X 


EH 


Eh 


Z 












z 




Z 


z 










a 


H 


Eh 


o 




< 


Z 


• 




o 


< 




z 




z 




z 






Eh 






z 


Z 










o 






u 




X 






CO 


Eh 






CO 


U 


o 


Z 


CJ 


u 


z 


£H 


Eh 




z 


z 




z 




Z 




O 


z 


H 










z 




O 


u 




% 


•ts* 




»o 


CO 


O 


w 


Z 


Z 


z 


t 


z 






Z 




t 


M 


z 


«» 


z 


Q 


O 


a 


Eh 


M 
















z 


X 






i> 




h3 


04 


z 


z 


W 


z 


z 


> 


H 


z 


z 


z 


O 


z 


O 


EH 


o 


z 


Z 


z 


z 




z 


O 










z 




<< 


z 




•-5 


<> 




-a: 


CO 


M 




Z 


z 




H 


to 


z 


z 


z 


E^ 


z 


z 




z 


z 


Eh 


1 


z 


O 




Z 










z 




£ 






M 


’Ll' 




u 


« 


O 




z 


Eh 


Q 




u 


z 


Eh 


z 




z 




z 


z 


z 


Z 


1 


CJ 


Z 




Eh 










z 




z 


M 




£H 










<< 


H 


z 




z 


Z 


z 


•a: 




z 


Z 




z 


o 


CO 


u 


z 


z 


z 


a 


• 












X 




a 


Z 




04 


-G* 




cn 


CO 


z 


z 


z 


O 


Eh 


O 


z 


> 


Z 


u 


z 


M 


z 


z 




u 


s 


z 


X 


a 


Eh 


Z 










£ 






z 




z 






H 


> 


w 


z 




Z 


«a 


H 


u 




H 




O 


z 


EH 




z 


z 


o 


o 


u 


< 


X 


X 










X 




z 


X 


Eh 


% 


'3> 






H 


> 


H 


z 


< 


z 


Eh 




rt: 




z 




*<3: 




o 


H 


z 


u 


H 






z 


X 










X 




o 


o 


X 




O' 




o 


CO 


c 


W 


o 




z 


CJ 


X 


z 


Z 


z 


CJ 


> 


z 


z 








EH 


z 


Q 


Eh 


X 
















z 


Z 


w 






z 


CO 




1 


a 


W 


u 


z 


z 


z 


Z 




z 




o 


z 


Z 


z 




< 


z 


z 


X 


u 










£ 




>• 


z 


z 


Q 






»< 


CO 


z 


X 


z 


z 




CO 


z 


H 


z 




M 


z 


z 




z 


o 


Z 


> 


z 


< 


o 


U 










O 




Eh 


z 


z 


O 


•G 






u 


o 


Eh 




O 


z 




> 


U 


u 


o 


O 


z 






z 




z 


Z 


z 
















Q 




M 


X 


£h 


z: 


-G* 




CO 


u 




Z 


z 




z 


z 


H 


Z 


u 


Eh 


O 


Eh 




a 


-a: 


z 


z 


z 




o 


z 












X 




Z 




z 








H 


ZD 


w 


o 




Z 


z 


M 


Z 


z 


z 




< 




z 


z 


> 


z 


z 


z 


z 


z 


z 


z 










< 




H 


z 


H 


M 


-Cl 




CO 


CO 


CO 


z 


w 


z 


Eh 




Z 


•pC 


z 


Z 




O 


z 


•< 




z 


z 


z 




E^ 


Eh 


z 










Z 




X 


o 






•& 




X 




o 




z 


M 




z 


z 






<c 


X 


H 


z 




< 


z 


u 


o 


X 


X 


C 


c 














< 




Z 


03 






h3 


CO 


z 


o 


o 


O 


CO 


z 


u 


a 




z 


z 




< 


z 


z 


z 


CJ 






z 


z 












z 




X 


z 


o 


CO 


'U> 




•< 


Eh 


z 


Eh 




Q 


M 


z 


u 


z 




o 




O 


> 


z 


z 


z 


ZD 


z 


z 


z 




a 






t « 




o 




o 


z 




Z 






z 


«; 


z 




o 


< 




H 


ZD 


z 


CJ 


z 


EH 


Z 




o 


H 




z 


o 


z 


o 


z 


z 


Eh 




Z 




z 




z 


X 


X 


% 






»< 


z 


z 


W 






Z 


Z 


Z 


H 


z 




u 


Z 


•< 


H 


u 


z 






z 


u 


o 


Eh 


X 




z 






z 


z 


£ 


z 


CD 








CO 




Z 






Z 


CJ 




z 


M 


z 


z 


CO 


z 




z 


z 


z 


z 


X 




z 


X 


z 




o 




>• 


X 




X 


M 


CO 


o 




CO 


z 


w 


O 


o 


Eh 


Z 


z 


X 


z 


Q 


o 


z 




z 


Z 


z 


Eh 


o 


z 


c 


z 




z 


Eh 




z 






z 


X 


z 


Eh 


n 






z 


CO 


z 




z 




z 


z 


z 


Q 


Z 




z 


Q 


M 


z 








z 


z 


M 


Eh 


s 


X 




z 




z 




z 




M 


% 






Eh 


o 


Eh 


ac 


w 


Z 


c 


o 






o 


Z 


z 


a 


u 


M 




Z 


z 


s 








o 


o 




o 




z 


z 


M 


Q 


Z 


a 












O 


CO 


H 


CO 




a 


z 


CJ 


c 


z 




z 


Z 


z 


M 


z 


z 


o 


Z 


X 


u 






z 






o 


z 


Z 


M 


CO 


•G* 




CO 


o 


z 


z 




W 




z 


z 


z 


z 


z 






z 


Eh 


X 


E- 


z 


z 


Q 




z 




z 












z 


Eh 


X 


w 


-& 




o 


EH 


o 


z 


Eh 


1 


o 


c 


CO 


EH 


z 


Eh 


z 


z 






Eh 


Z 


s 




ZD 


z 




z 


z 




z 






EH 


£ 


U 


••a; 


cn 


■» 








z 




z 


X 


a 




X 








z 


z 




H 




z 


z 


z 


z 


o 


z 


M 


Eh 




c 




s 




rt: 


Z 


X 


M 


'S' 




Eh 


z 




CO 


z 


£H 


ZD 


z 


z 


Eh 


Q 


Z 


<< 


z 


z 


Q 


Z 




z 


X 


z 




z 




< 








£ 


O 


X 


z 










Z 


c 


03 


z 


M 


Z 


z 


z 


!< 


Z 


Z 


Z 






z 




Z 


X 




Eh 


z 


z 


z 


Z 


Z 




z 




X 


z 


z 


X 


Z 


CO 






-a: 


z 


W 


w 


Z 


n 


CO 


o 


z 


Z 


c 


Z 


z 


z 


CJ 


Z 


z 


z 


Q 






PH 


-c 


O 






z 




Q 


z 


z 


z 


z 








CO 


o 


lO 


z 


1 


z 


z 


M 


<< 


M 




Z 


z 


M 


z 


-< 


X 


z 


z 


z 


z 


•< 


M 


z 


Z 




c 














04 


42* 




z 


o 


z 


w 


X 






Eh 




*< 


> 




z 


rt: 


Eh 


M 


z 


o 


E^ 


z 


z 


z: 


Z 


z 


o 








1 


1 


1 


1 


1 


z 








z 


z 


Eh 


Eh 


X 


z 


C 


z 


o 


z 


Z 


z 


z 


Z 


S 






o 


s 


M 


H 


< 


z 


z 




Q 




1 


1 


1 


< 


1 


#tc 






CO 


z 




z 


Z 


z 


z 


> 


H 


c 


H 


M 


z 




H 


o 


z 


Q 


z 


H 


EH 


Eh 


> 








Z 














cn 


i;J 




z 




CO 


H 


O 




Eh 


z 






X 




u 


z 




z 


z 


Z 


z 


£H 


z 


Z 






z 




z 




















Eh 


z 






Z 


z 




z 


Z 


z 




z 


u 


z 


Z 


M 


o 


a 


X 




z 


z 


z 


z 


z 




X 














CO 








M 


o 


Z 




z 




z 


z 


o 


X 


o 


ZD 


H 


O 


Z 


H 


z 


z 


Z 






o 


z 


s: 


















z 






CO 




a 


O 


Z 


M 


• 


z 


z 


H 


fJO 


H 


z 


z 






Eh 


O 




Z 


z 


z 


£H 


< 






z 














M 


«• 




M 


z 


z 


a 


Z 


X 


CO 


o 


z 


EH 




EH 




z 


Z 


Z 


Z 


u 


z 


H 


X 


z 


c 


X 


z 




z 








X 






Eh 


<{3- 








z 


z 


M 


hO 


z 




< 


<a: 


z 


< 






z 


> 


z 


z 


z 


EH 


Eh 


X 


z 


Oi 


z 




z 








o 


X 




3 






CO 


w 


CO 




Q 


CU 


z 


Q 


z 


> 


z 


> 


f 


z 




M 


H 


z 


Eh 


Z 




a 


X 


z 


X 




X 








z 


z 




o 


* 




hJ 


z 


z 


z 


Z 


H 


o 


z 




z 


H 


Z 


z 


z 


CJ 


Eh 


Z 






Z 


Z 


z 


z: 




•a: 










z 


z 


X 


X 


03 






04 


Eh 






*< 


Eh 


z 




z 


z 


Z 


Z 


ZD 


u 




C 


Eh 


z 


Z 


U 


Z 


u 


X 


z 


H 




H 




••a: 


< 


X 


z 


z 


03 






s 




w 


w 


Z 






Z 


CJ 


z 


Z 


z 


Z 


CJ 


z 


C5 


Z 


z 


•< 


Z 


Eh 


o 


u 




z 




z 














ZD 






c 


X 


z 


z 


X 


Z 


Z 


z 


< 


m 


-< 


z 


••a: 


ZD 


z 


z 


M 


z 


X 


z 


Z 


z 


u 


z 


c 


















CO 


<f 




CO 


Z 


Eh 


EH 


w 


z 


H 


Q 


z 


o 


z 


o 


> 


z 


Eh 


z 


Q 


z 


Eh 


z 




z 




£ 


> 




> 
















V* 

i* 

i? 


•& 


•G* 


<3* 


i? 




42 


42 


42 


42 


«• 


Ci- 


•a.* 


42 


42 




42 


42 


<h 


o 


<1 


<i> 


42 


42 


42 


42 


42 


42 




42 


42 


42 


42 


42 


42 




<5* 


•wi' 










42 


c; 


42 




42 


M 




i2 




42 


<2 


w 


42 


w 


42 


42 


'C2 


•e2 


42 




'i2 










42 




•S' 


42 




U 


U 


U 


U 


u 


u 


U 


u 


U 


U 


U 


u 


U 


U 


u 


U 


U 


U 


U 


u 


U 


U 


U 


CJ 


U 


U 


U 


U 


u 


U 


u 


U 


u 


u 


U 



115 







•{n» 


4H* 












42. 


<£;■ 








<n. 




«£> 






<2- 


























,J2* 


<3* 














•2* 
































































































o 






























































































































CO 


























■2 






















H 










W 












z 














































X 


CO 










X 




































U 










04 


z 










> 














•2 






















M 










Z 


o 










H 




































EH 










i< 












EH 






X 








O 






















C/5 








CO 


CO 


z 










•< 








z 






O 




% 


















O 








IX 




o 










> 






M 


X 










X 




EH 














Z 




Q 




h4 


04 


z 










H 






Z 


> 






§ 




X 




Z 














o 




Z 




M 


l< 


X 










Z 




z 


EH 


M 










M 




M 














c 




03 




EH 


z 












X 




z 




CD 










EH 




Z 














M 






CO 


Z 


w 


z 










X 




o 


X 












Z 


% 


Z 














Q 








w 




o 




Z 










H 


u 


< 






■J* 




Z 


X 


M 






•k 












03 


w 


u 


04 


H 




X 


X 




X 




£h 


c 






























Pu. 




O 


Q 


z 


z 


£h 




X 


X 




z 


z 


-a: 


X 


z 




















X 








O 




EH 


O 


w 


< 






M 


z 




< 


X 


> 




o 




















X 












C 


Z 


04 


CO 


> 




EH 


z 


Eh 




EH 


z 


z 


X 






■to- 














M 








EH 




03 








03 




z 


c 


Z 


z 


< 


X 


X 






















EH 








r> 




W 


z 


2^ 


z 


w 




X 


z 


X 


z 


z 


z 


z 


z 




















Z 








04 




Z 


CD 


CD 


CD 


CO 




u 




EH 


o 


H 


X 


z 


o 




















Z 


X 






EH 




W 


z 


z 


z 


z 




z 


z 




H 


£H 


o 


X 


M 






















X 






3 




O 


O 


O 


o 


o 


£h 


X 


H 


O 


EH 


Z 




u 


EH 






'S' 














X 


H 






O 






03 


z 


03 




z 


z 




EH 


< 


X 


X 


u 


-s5 




















X 


EH 










03 


z 


z 


z 


z 


X 




Z 




:> 




z 


z 


> 










X 










X 


Z ^ 






« 




W 


EH 


EH 


Eh 


o 


Eh 


X 


z 


Z 


z 


X 


c 


z 


z 






2 




X 




EH 






o 


Z X 






o 


CO 


CO 








Q 




o 


o 


X 


X 


X 


z 




X 


X 




‘2’ 




X 




Z 






z 


^ X 


A 






H 


s 


CD 


CD 


CD 


Z 


O 




H 


X 


z 


M 


X 


H 


z 


X 




•2> 




o 


% 


Z 






%X M 


o 






Q 


SD 


Z 


Z 


Z 


c 


EH 


z 


£h 


H 


X 


Eh 


X 


X 


X 


X 




-»r*. 

V* 




z 




H 






rn 


X Eh 


r* 




cn 




z 


H 


H 


M 


« 




o 




H 


o 


Z 




X 


o 


o 
















X 




X Z 


w 




w 


CC 




U 


CJ 


U 




Z 


z 


> 


Z 




X 


X 


Z- 




z 




<> 












o 


EH 


O Z 


Z 




£h W 


z 


Z 


z 


z 


EH 


X 


z 


Z 


X 


X 


u 


o 


o 


X 






<> 












z 


< 


z ^ 


z 




w 


£H 


o 


w 


w 


w 


CJ 


X 


X 


X 


a 


z 


z 




z 


z 


X 
















z 


X 


X 


X 


o 


c 


Z 


a 


z 


z 


z 


w 


Q 




z 


z 


•ec 


X 


X 


z 


*< 


z 
















X 


z 


X X 


X 


o 




M 


z 


o 


o 


ot 


X 


O 


z 


X 


X 


z 


z 


X 


X 


z 


<3 
























w 


03 


CX 




w 


IX 


X 


X 


Z 


X 


o 


z 


X 




z 


X 


X 


























EH 


C 


o# 


03 


CO 


CO 


CO 


z 




z 






X 


X 


z 




X 


X 
























< 


04 














X 


<3 


X 


X 




X 


< 


X 




X 
























U 




x 


03 


X 


X 


X 


O 


o 


z 


o 


c 


X 


M 


z 


o 


X 


w 




























o 


O 


w 


w 


w 


£H 




o» 






o 


EH 






X 


Eh 




-3* 
















«b 




Q 


o 


IX 


X 


Q 


Q 


Q 




z 


z 


z 


z 




< 


o 


z 


EH 


z 












% 












z 


(X 






z 


Z 


z 


X 


X 




X 


X 


X 


X 


X 


X 


u 


X 








z 




X 








X 






£H 


z 


Q 


M 


H 


H 


X 


z 


z 


X 


X 


<3 


z 


z 


X 


X 


CJ 








X 


% 


z 








X 




s 


z 


03 


W 








Q 


z 


<3 


z 


z 


z 


z 


X 


z 


z 


z 




■e> 




z 


n 


M 








M 




w 


O 


CO 


M 


O 


O 


o 


Z 


z 


X 


z 


z 


z 


z 


z 


z 


X 


X 




















EH 




Q 


o 


Q 


CO 


Q 


Q 


o 


M 


£ 


Z 


z 


z 


*< 


o 


z 


z 


X 


z 




















z 


























































z 




1 


f 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


























1 


\ 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


f 


1 


1 






















X ^ 


Nv 
























































X ^ 


z 






















































% 


X ro 


X 








































■& 
















o % 


X 






















































X 


z z 


M 






EH 


















X 








z 










43* 


% 




Z 








^ CQ 


% X 






Z 


Q 










X 






X 








X 








* 


CJ' 


X 




X 








m X 






EH 


M 


W 










X 






M 




EH 


z 


u 










z 


X 




u 






z 




X 




« 


03 


IX 








Q 


Q 


X 


z 


Eh 


z 


<5 


z 


u 








*2* 


X 


X 




u 


43* 




z 


X z 


ZT Z 




04 


04 


CO 








Z 


O 


z 


X 


Z 


X 


X 




z 








Cl- 


CD 


z 




X 


& 


z 


•< 


Z X 




H 


H 


H 


M 






X 


H 


z 


z 


z 


Z 


o 


z 


z 


z 


X 






'S' 


X 


M 


H 


z 


X 


<3 


z 


z o 


< z 








































C 


Eh 








-<3 








z 










































Z 








X 








o 










































H 








z 








u 








































C 

•S' 






42- 


4S- 






Cl 


43- 


43 


O 




-S* 






C 






<& 


e 




•2 






"TJ 

>r 
































«• 






*1? 


il’ 


# 


V? 








*2 












<4* 




•2 




















U 


u u 


U 


U 


u 


U 


U 


u 


u 


u 


U 


u 


u 


U 


u 


u 


u 


u 


u 





















116 



IPRT, ICONIN, ICNOUT, IPR INT 

IDSK (30) 

AN = FLOAT (NBR) 

DO 10 K=1,NPTILE 

EXPS (K) = FX(K) AN 



GENERATE PSEUDO SAMPLES 
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IN THE CASE OF INTEGER VALUES OF EXPECTED 
DEMAND, SUCCESSES ARE ACCUMULATED UNTIL 
THE EXPECTED NUMBER OF SUCCESSES IS 
REACHED. 






I 



IF (FLOAT (SUCCES) .GE. EXPS(K)) GOTO 40 
SUCCES = SUCCES + 1 
CONTINUE 
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SUBROUTINE OUTPUT (MODEL, NPT ILE, NREPS , FX , IC , PHAT , MSE , CAT) 
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WRITE(IPRINT,610) DISTR (J) , (FX (K) , K=1,NPTILE) 

TMSE =0.0 
DO 10 K=1,NPTILE 

MU(K) = PHAT(I,J,K) / TRIALS 

STDV (K) = SQRT(MU(K) « (1.0 - MU(K))/ TRIALS) 
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